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The  kinetics  and  mechanisms  of  interactions  of  nitrogen  and  carbon 
monoxide  with  liquid  niobium  were  investigated  in  the  temperature  range  of  2700  to 
3000  K in  samples  levitated  in  NJAr  and  CO/Ar  streams. 

with  the  standard  enthalpy  and  entropy  of  solution  of  - 236.4  ± 23.3  kJ/moI  and  - 
65.3  + 8.3  J/K*mol,  respectively.  Accordingly,  the  nitrogen  solubility  (C*.  in  wt%) 
is  expressed  as 


where  T is  the  temperature  in  K,  and  Pw  is  the  nitrogen  partial  pressure  in  atm. 

on  and  desorption  processes  were  found  to  be  second- 


adsorbed  nitrogen  molecules  inio  adsorbed  atoms.  The  rare  equations  for  nitrogen 
absorption  and  desorption  are  given  by 

( -T?  >_  “ 2 M°-10’  “P  <-  ■ 3 T°°  ) -Jr  «ft*  ‘ <V) 

( -t? >„  = '•6*°*10’ «p<-  61t°  > 4- (ft.*- <V> 

The  carbon  and  oxygen  dissolution  in  liquid  niobium  (tom  CO  gas  is  also  an 
exothermic  process  and  the  solubilities  of  carbon  and  oxygen  (C^,  Co,  in  at%  ) are 
related  to  the  temperature  and  the  partial  pressure  of  CO  (P„  in  atm  ) as: 

In  CQ  + In  Co,  - In  P„  = 2.32-10YT  - 2.35 
The  reaction  CO  ->  [Q  + [01  along  with  the  evaporation  of  niobium  oxide 
takes  place  during  C and  O dissolution,  whereas  C and  O desorption  occurs  via  CO 
evolution  only,  especially  when  excess  carbon  exists  in  the  liquid.  Absorption 
/desorption  of  carbon  and  oxygen  in  liquid  niobium  are  second  order  processes  with 
the  first  order  dependence  on  C and  O concentrations  (Cc,  Co  in  al%)  respectively, 
implying  that  the  overall  reaction  rate  is  controlled  by  the  substep  of  dissociation  of 
adsorbed  CO  moleculcsfor  recombination  of  C and  0 in  the  case  of  desorption). 
The  rate  equations  for  C and  O dissolution  and  desorption  are  found  to  be 

(dftdlU  = (l/0.76)(dCo/dt)„  = 12.1(A/V)exp(-26700/T)(Cc,Cu,  - QXy 


(dQ/dt)»  = (dCo/dt)*,  = 3.80*lO‘(A/V)cxp(.47500/T)(CClCo,  - Q.  CJ 


CHAPTER  I 
INTRODUCTION 

A unique  combination  of  properties  such  as  relatively  low  density,  high 
melting  point,  high  superconducting  transition  temperature,  low  thermal  neutron 
absorption  cross  section,  and  good  fabricability  makes  niobium  and  its  alloys 
attractive  candidate  materials  in  a variety  of  applications.  Its  use  ranges  from 
superconductors,  to  high  temperature  structural  applications,  and  to  core  materials  for 
nuclear  fusion  reactors.  However,  because  of  the  great  affinity  of  Nb  and  niobium 
base  alloys  for  interstitial  elements  like  H,  N,  O,  and  C,  these  materials  react 
readily  with  the  reactive  gases  containing  the  interstitials,  causing  striking  changes  in 
their  physical  and  mechanical  properties! 53Ang,  63Des,  66Wis,  70Els,  71Mey, 
72Szk,  73Kle,  73Gio,  81Eng).  For  example.  1 at.%  nitrogen!"  0.15  wl%)  dissolved 
in  niobium  increases  its  hairiness  by  about  270%  and  its  electrical  resistivity  by 
about  30%(Fig.  l.l)[63Gcb,  66Geb].  The  effects  of  interstitials  on  various  properties 
of  pure  Nb  are  summarized  in  Table  1.  Therefore,  a strict  control  of  the 
interstitial  content  is  important  for  the  reliable  use  of  these  materials  in  advanced 
technological  applications. 

Conventional  production  route  of  high  purity  Nb  consists  of  (1)  the 
extraction  of  NbjOj  from  niobium  ores  through  leaching  or  chlorination/distiiluiion, 
(2)  carbothermic  or  aluminothermic  reduction  of  Nb,0,,  (3)  consolidation  process  to 
form  the  electrodes  through  sintering  or  remelting  of  Nb  powder,  and  finally  (4) 


Fig.  1.1  Effect  of  dissolved  gases  on  the  mechanical  and  elec 
of  recrystallized  niobiomI63Geb.  66Geb]. 


ctrical  propertie 


Table  1.  Properly  changes  of  Nb  due  I 


dissolved  H,  N,  O,  and  C at 


specific  property  Nb  - H 

changes 


Nb  - O Nb  - C 


Lattice  parameter  2.3x10*  6.0* 104  6.1x10* 

Aa/Ac,  nnt/at.% 


supplementary  electrical  0.66  4.1  4.1 

resistance,  Ap/Ac 

(iflcm/at.% 

Magnetic  susceptibility  - -5.3*10*  -5.7*10* 

decrease  AyJAc, 

cm7g-at.% 

Hardness  increase,  6.9  ISO  120 

Vickers  hardness  number 
per  1 at.% 


Hardening  coefficient  - 591  382 

Ao/Ac,  MN/m’  al.% 


4.3 


after  ref.|80Fro,  SISchl  | 


refining  through  electron  beam  remelting  in  vacuum  (81Sch,  85Asf,  85Sch].  Among 
the  above  4 stages,  purification  of  Nb  is  mostly  achieved  during  electron  beam 
remelting  in  vacuum  (EBM)  and,  in  part,  during  consolidation  step  in  high  vacuum 
through  degassing  of  nonmetallic  elements  and  volatilization  of  metallic  impurities. 

In  EBM,  the  electron  beam,  accelerated  in  a high  potential,  is  focussed  on  to 
the  Nb  ingot  to  be  melted,  which  drips  down  into  a water  cooled  mold.  The 
solidified  ingot  is  then  withdrawn  downward  on  a base  of  the  same  material.  Since 
this  process  is  performed  in  high  vacuum,  non  metallic  impurities  such  as  H,  N,  O. 
and  C,  and  volatile  metallic  impurity  like  Al  are  removed  to  a great  extent  via 
degassing  and  vaporization  process.  Table  2 shows  the  typical  impurity 
concentrations  before  and  after  EBM.  In  some  cases,  a special  purification  stage 
such  as  electrolytic  refining  in  molten  salt  can  be  incorporated  before  EBM  in  older 
to  eliminate  the  less  volatile  elements  like  Ta  and  W.  On  the  other  hand,  Nb  in 
solid  state  can  be  purified  in  ultra  high  vacuum  at  high  temperature  above  2000  °C, 
but  purification  kinetics  is  much  slower  than  that  in  liquid  state  because  transport 
from  the  interior  to  surface  in  solid  state  obeys  the  law  of  diffusion  while  in  liquid 
metal  transport  of  impurities  take  place  via  hydrodynamic  flow.  For  example,  it  is 
repotted  by  K.  Schultz  and  cl  al.  181Schl.  81Sch2J  that  an  annealing  time 
longer  than  one  week  is  required  to  reach  the  ppm  level  for  a niobium  rod  of 
10  mm  diameter  at  2400  K under  the  nitrogen  partial  pressure  of  5*10"  mbar. 
Thus,  EBM  has  been  used  widely  in  the  industries  as  a major  purification 
step  in  the  production  route  for  high  purity  niobium  since  1950s. 


Table  2. 


atmosphere  : 3*10*  mbar 

ingot  melting  spccd(kg/h):  8.25(lst  run),  12.85(2nd  run), 
59  (3td  run) 


Element  ATR** ***  ingot  1st  run 


1450  1450 


* after  ref.[81Schl] 

*'  ATR  : Aluminothermic  reduced 

***  ND  : None  detected  («  10  ppm) 


The  first  fundamental  results  concerning  Nb  and  gas  reacdons  ai  high 
temperature  and  low  pressures  were  reported  around  1965I65Froll-  Since 
then,  a considerable  amount  of  research  has  been  earned  out  to  understand 
the  interaction  behavior  of  interstitials  in  Nb  but,  in  most  cases,  limited  to  the 
solid  state.  The  high  melting  point  of  Nb  (2468  °Q  precludes  the  use  of  the 
conventional  experimental  techniques  used  to  study  such  processes.  As  a 
result,  there  is,  at  present,  quite  a lack  of  fundamental  understanding  of  the 
interaction  processes  occuring  between  the  gas  phase  and  liquid  Nb.  which  is 
essential  for  efficient  refining  and  control  of  the  interstitial  content  in  the  melt. 

The  aim  of  this  research  is  to  attain  a comprehensive  understanding  of 
mechanisms  and  kinetics  of  absorption  and  desorption  behaviors  of  N,  C,  O 
in  liquid  Nb.  The  emphasis  of  the  study  is  to  determine  the  solubilities  of 
each  interstitial  in  liquid  Nb,  reaction  mechanisms,  rate  controlling  step,  and 


The  limitations  of  conventional  techniques  were  circumvented  in  this 
study  by  carrying  out  the  experiments  in  an  electromagnetic  (EM)  levitation 
apparatus  [65Fro3,  65Pei,  84Eth].  Electromagnetic  levitation  is  particularly 
suited  for  the  present  study  since  this  containerless  processing  method 
minimizes  impurity  pick-up  from  extraneous  sources.  Also,  the  high  melting 


the  advantages  of  rapid  heating  under  various  gas  atmospheres  followed  by 
rapid  quenching  if  desired.  The  sample  can  be  cycled  through  desired 
e ranges  for  long  periods  of  time,  thus  offering  a way  to  monitor 


process  under  various  processing  conditions. 


CHAPTER  2 

REVIEW  OF  LITERATURE 


In  reecni  commercial  Nb  production  lines,  purification  of  niobium  is  mostly 
conducted  by  EBM  in  high  vacuum.  Since  the  reduction  of  interstitial  impurities  in 
this  stage  is  achieved  through  the  gas-metal  reaction  at  high  temperature  and  low 
gas  pressure,  thermodynamics  and  kinetics  studies  associated  with  it  arc  both 
important  to  understand  the  purification  process.  Thetmodynamics  gives  the 
informanon  not  only  of  the  direction  of  reaction  but  also  of  the  final  composition 
attainable  under  a given  processing  condition.  Kinetic  data  will  provide  the  idea  of 
how  fast  the  reaction  can  reach  the  equilibrium  and  how  reaction  rate  can  be 
adjusted  toward  the  desired  direction. 

Therefore,  in  the  sections  to  follow,  the  thermodynamics  associated  with  gas- 
metal  reactions  will  be  discussed  first,  followed  by  the  discussion  of  mechanism  and 
kinetics  involved. 


Thermodynamics  of  Solution 

When  a metal  reacts  with  gases  such  oxygen,  nitrogen,  and  CO,  it  depends 
upon  temperature  and  partial  pressure  of  the  gas  whether  the  reaction  leads  to  the 
compound  formation  such  as  oxides,  nitrides,  and  carbides,  or  to  the  interstitial 
solution  in  metal.  Below  the  dissociation  pressure  of  the  compound  at  a given 
temperature,  the  gas  species  are  dissolved  in  metal  in  the  fotm  of  interstitial 
solution,  whereas  beyond  that  pressure  compounds  start  to  form  on  the  surface  of 
metal.  Discussion  will  be  stoned  with  solution  reaction. 

Suppose  that  a pure  metal  is  in  equilibrium  with  its  vapor  as  follows. 


(2.1) 


The  free  energies  of  Ihe  Iwo  phases  in  equilibrium  are 

G„  = G_  (2.2) 
where  G„  and  G„  are  the  free  energies  of  a melal  and  its  vapor,  respectively.  For  a 
constant  temperature,  the  change  in  free  energy  of  vapor  phase  due  to  its  pressure 
change  is  calculated  using  the  following  equation. 

dG,  = VdP,  (2.3) 
where  V is  the  volume  of  the  melal  vapor  and  P is  the  vapor  pressure.  Integration 
of  equadon(2.3)  gives 

G„  = G„°  + RT  In  0VP.°)  (2.4) 
where  P,  is  the  equilibrium  vapor  pressure,  R is  the  gas  constant,  T is  the 
temperature,  and  ° indicates  the  properties  at  standard  pressure.  If  the  pressure  is 
given  in  atm  and  the  standard  stale  is  taken  l atmosphere,  the  (fee  energy  of  vapor 
pressure  is 

G„  = G„°  + RT  In  P.  (2.5) 
On  the  other  hand,  for  the  metal,  the  tree  energy  is  given  by 

G»  = G°„  (2.6) 
because  the  change  in  free  energy  caused  by  the  change  in  pressure  is  negligible. 
According  to  equation(2.2).  it  can  be  seen  that 

G„  = G°„  = G„  = G0„  + RT  In  P. . (2.7) 
The  change  in  the  free  energy  due  to  vaporization  is  therefore 

G°„  - G°„  = AG°  = - RT  in  P.  . (2.8) 
This  equation  shows  how  the  standard  free  energy  change  for  vaporization  is  related 


the  equilibrium  vapor  pressure. 


So  far  only  equilibrium  belwccn 


considered.  However,  in  many  meiallurgical  processes,  meiallic  solutions  of  several 
components  are  usually  involved.  Thermodynamic  treatment  of  these  solutions  is 
dealt  with  the  partial  free  energy,  which  is  defined  as  the  change  in  the  free  energy 
of  the  solution  when  one  mole  of  a substance  is  dissolved  in  a very  large  amount 
of  solvent.  In  discussion  regarding  the  equilibrium  between  pure  metal  and  its  vapor 
pressure,  it  was  shown  how  the  free  energy  of  the  gas  or  vapor  phase  is  varied 
depending  upon  the  change  of  pressure  from  P,°  to  P..  In  a similar  manner,  the 
change  in  the  free  energy  of  solution  when  one  mole  of  substance  i is  added  can  be 
derived  by  substitution  of  the  partial  pressure  P,  and  P,°  for  P,  and  P,°,  respectively, 
leading  to  the  following  expression. 

AO,  = G,  - G,°  = RT  in  (P/P,°)  = RT  In  a,  (2.9) 

where  G,  and  P,  are  the  free  energy  and  partial  pressure  of  substance  i,  respectively, 
and  G,°  and  P,°  indicate  the  same  properties  at  standard  state,  respectively.  The  ratio 
P/P,"  is  called  the  activity  of  substance  i,  a,,  in  the  solution.  If  X,  is  the  molar 
fraction  of  the  component  i in  the  solution,  the  following  relation  exists  between  its 
activity  and  the  concentration  X,: 

a.  = Y.  X,  (2.10) 

where  y is  the  activity  coefficient  of  clement  i.  Equation(2.9)  becomes  therefore 

AG,  = RT  In  y + RT  In  X, . (2.1 1) 

For  an  ideal  solution,  y = 1,  a,  = X,.  Such  a solution  has  no  interactions  of 
the  components,  obeying  Rault's  law.  If  y > 1,  there  is  a tendency  for  separation 
between  solute  and  solvent  atoms  and  a reduced  solubility  is  observed  for  a given 


partial  pressure  of  component  i.  Finally,  if  y, 


compound  formation  between  the  solvent  and  dissolved  substance,  and  a higher 
solubility  is  found  than  an  ideal  solution.  Regarding  the  temperature  dependence  of 
the  activity  coefficient,  temperature  increase  of  a non-ideal  solution,  in  principle, 
results  in  a decrease  in  the  extent  to  which  the  activity  coefficient  deviates  from 
ideal  solution.  In  other  words,  if  y,  > 1,  then  an  increase  in  temperature  decreases 
the  value  of  Y,  toward  unity,  increasing  the  solubility  at  a given  partial  pressure, 
whereas  if  Y.  < 1.  then  the  value  of  y,  increases  toward  unity  with  increasing 
temperature,  resulting  in  the  reduced  solubility.  It  can  be  therefore  predicted  that  if 
y is  greater  than  unity,  the  solution  reaction  is  endothetmic,  while  if  y,  is  less  than 
unity,  it  is  exothermic. 

Now,  turning  the  attention  to  a dilute  solution,  it  should  be  noted  that  there 
always  exists  a composition  range  with  a constant  activity  coefficient,  y.  If  the 
constant  value  of  the  activity  coefficient  is  set  to  be  y*,  the  partial  pressure  of 

P,  = P“,  a,  = - k,,  X,  (2.12) 

which  is  the  Henry’s  law  and  k,,  is  Henry’s  constant.  Fig.  2.1  shows  a typical 
diagram  of  the  partial  pressure  of  the  components  of  a binary  alloy  system  with 
positive  deviation.  When  Rault's  law  is  obeyed,  the  partial  pressure  is  linearly 
proportional  to  molar  fraction,  lying  on  the  straight  lines.  Henry’s  lines  are  the 
tangents  to  the  activity  curve  for  XA  or  X,  ->  0.  respectively.  As  may  be  seen  from 
the  diagram,  yA  and  Y„  are  constant  only  over  a very  small  range  of  concentration. 


The  partial  pressure  curves  also  demonstrate 


Examples  of  (he  relationship  between  vapor  pressure  and  concentration 
in  a binary  system  A-B  with  a positive  deviation  from  Rault's  Law 
|71Win2). 


solubility  is 


Henry's  activity  based 


Ihc  pure  subsiance  as  "standard  state."  N 
Henry's  line,  a”,  is  used,  which  is  defined 
■ f,  C, 


where  Henrian  activity  coefficient,  f„  becomes  unity  when  the  concentration  of  the 
specie  i,  Q,  ->  0.  The  concentration  C,  is  conventionally  expressed  in  weight  %. 
Henry’s  activity,  aH,  is  therefore  equal  to  the  concentration  of  the  dissolved 
substance  C,  in  weight  % in  ideal  dilute  solutions. 


Since  most  of  gases  are  dissolved  in  metals  almost  solely  in  their 
state,  the  solution  reaction,  for  instance,  for  nitrogen  can  be  represented  as 
1/2  (NJ  .[N)„ 
with  the  equilibrium  constant. 


K = — (Pv.)1" = exp(  -AG7RT)  (2.15) 

where  a»  is  the  activity  of  nitrogen,  P*,  is  the  partial  pressure  of  N„  f is  the 
Henrian  activity  coefficient,  C.-  is  the  nitrogen  concentration,  AG°  is  the  free  energy 
change  for  solution  reaction  of  one  mole  of  nitrogen  atom,  R is  the  gas  constant, 
and  T is  the  temperature.  For  low  concentrations  where  the  activity  coefficient  is 
almost  constant,  solubility,  C„,  is  proportional  to  the  square  root  of  partial  pressure, 
V.  for  diatomic  gases.  This  is  the  well  known  Sicvert's  law  and  can  be  expressed 


C»  = k(T)  (P* 


(2.16) 


where  k(T)  is  iemperalure  dependent  Sievcns  constant.  From  the  equation 
following  relationship  can  be  derived  as 


(2.17) 


C,  = (1/0  (P»)'° 


where  AH0  and  AS°  are  solution  enthalpy  and  entropy  change,  respectively.  Hence,  f 
is  unity  in  dilute  solutions  where  Henry's  law  is  valid  and  standatd  state  is  taken  as 
1 wl%  solution.  Therefore,  comparing  cquation(2.16)  with  (2.17),  it  can  be  seen 


Substituting  the  equation  (2.18b)  into  equation  (2.16),  it  is  found  that  solubility,  C„ 
is  related  to  the  partial  pressure  and  temperature  as 


Equation  (2.19)  indicates  that  the  solubility  at  a constant  pressure  is  increased  with 
rising  temperature  if  AH°  is  positive,  i.e„  if  the  solution  reaction  is  an  endothermic 
process,  while  the  solubility  is  decreased  with  the  increase  of  temperature  if  AH°  is 
negative,  i.e.,  if  the  solution  reaction  is  exothermic.  On  the  other  hand,  the  solubility 
at  a given  temperature  is  increased  with  increasing  partial  pressure. 

In  practical  application  of  solution  reaction  of  the  intetstitials  in  metal  such 
as  vacuum  degassing,  it  is  necessary  to  examine  how  the  solubility  varies  with  the 


(2.18b) 


(2.19) 


change  of  partial  pressure  and  temperature.  Fig.  2.2  shows  the  pressure  dependence 
of  the  solubility  of  various  gases  (also  those  of  sulphur  and  phosphorous)  in  the 
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Fig.  2.2  Solubilities  of  gases  in  liquid  iron  at  1600  "C.  in  liquid  nickel  at 
1450  °C  and  liquid  chromium  at  1745  °C.|71Win2] 


liquid  iron,  nickel,  and  chromium  |71Win2|-  II  may  be  seen  dial  only  for  hydrogen 
and  nitrogen,  the  equilibrium  pressures  are  high  enough  even  at  low  concentrations 
that  a vacuum  removal  of  the  elementary  gases  is  feasible  by  lowering  the  pressure. 
This  fact  is  quite  important  in  vacuum  degassing  because  degassing  is  possible 
without  a significant  loss  of  metal  only  when  sufficiently  high  molar  concentrations 
of  the  gas  exists  in  the  vapors  coming  out  from  the  metal. 

Fig.  2.3  shows  the  temperature  dependence  of  equilibrium  pressures  of 
oxygen  at  a given  oxygen  level  of  100  ppm  for  iron,  nickel,  and  copper.  The  slope 
of  these  curve  is  markedly  smaller  than  that  of  the  vapor  pressure  curves  of  the 
solvents.  It  may  therefore  be  predicted  that  no  improved  degassing  conditions  can 
be  obtained  for  oxygen  by  raising  the  temperature. 

Temperature  dependence  of  the  solubility  of  gases  varies  greatly  from  metal 
to  metal.  Fig.  2.4  shows  the  nitrogen  solubilities  in  group  Vb(Tn,  Nb)  and  VIb(W, 
Mo,  Cr)  solid  metals.  Solid  lines  indicate  the  solubility  lines  whereas  dot-dashed 
lines  represent  the  phase  boundary  between  a solid  solution  and  compounds.  With 
increasing  the  temperature,  nitrogen  solubilities  increase  for  VIb  metals  but  decrease 
for  Vb  metals.  This  indicates  that  the  sotudon  reaction  for  Nb  and  Ta  is 
exothermic! AH  negative)  whereas,  for  group  VIb  metals,  it  is  endothermic! AH 
positive).  At  temperatures  between  1600  and  2400  °C,  niobium  and  tantalum  show 
rather  high  nitrogen  solubilitics(bctween  0.1  and  10  at9f>)  at  relatively  low  nitrogen 
pressure  of  104  to  1 tore.  The  terminal  solubililyfdoi-dashcd  lines)  of  several  at% 
can  easily  be  reached  at  temperatures  of  about  2000  °C  and  below.  In  contrast,  the 
nitrogen  solubility  in  Mo,  W,  or  Re  is  very  low(10’  at%)  even  at  pressures  of  about 


[ 

{ 

r 

/ 

{ 

t 

/ 

/ 

/ 

■7 

1 

_L 

looo  isoo  1400  1600  teoo  2000 


Temperature  dependence  of  the  oxygen  pressure  of  copper,  iron,  and 
nickel  containing  100  ppm  oxygen  and  of  the  vapor  pressure  of  the 
pure  meials.[71Win2) 


Fig.  2,3 


Solubilities  of  nitrogen  in  solid  Nb,  Ta,  Cr, 


and  Re.[72Frol] 


1 aim  N,.  Based  on  those  facts,  it  can  be  stated  that  degassing  treatments  in 
vacuum  are  most  effective  for  group  Vb  metals  at  very  high  temperatures,  whereas 
for  group  VIb  metals  the  annealing  temperatures  should  be  low  as  far  as  kinede 
consideration  allows. 

For  oxygen,  solid  solution,  as  a rule,  is  expected  only  in  very  low 
concentration  ranges  because  of  higher  affinity  to  metal  compared  with  that  of 
hydrogen  or  nitrogen.  In  such  a limited  composition  range  the  deviation  from 
Henry’s  line  is  usually  very  small,  as  demonstrated  in  Fig.  2.5,  in  which  solubility 
of  oxygen  in  iron,  nickel,  and  their  alloys  is  plotted.  In  addition,  the  maximum 
attainable  equilibrium  pressure  for  a given  temperature  corresponds  to  the 
decomposition  pressure  of  the  oxide,  which  is  in  equilibrium  with  the  melt  at 
saturation.  Therefore,  the  terminal  solubility  of  gas  atom  in  the  solid  or  in  the  liquid 
can  be  determined  by  thermochemical  stability  of  oxides  with  respect  to  the  solution 
phase.  Consider  the  reaction  of  oxide  formation  for  the  metal-oxygen  reaction  as 

1/2  0,(g)  + mM  -a  M.0  AG,0  (2.20) 

and  the  solution  reaction  as 

1/2  0,(g)  -*  0(in  M)  AG,°  . (2.21) 

The  equation  for  the  oxide  formation  from  saturated  solution  is  derived  by 
subtracting  equation  (2.20)  from  equation  (2.21). 

M.0  ->  0(in  M)  + mM  AG,°  (2.22a) 

AG,0  = AG,”  - AG,° . (2.22b) 

Equilibrium  constant  for  reaction  (2.22a),  Kp,  is  given  by 


Fig.  2.5  Oxygen  pressure  as  a function  of  the  oxygen  conieni  of  pure 
pure  nickel  and  several  iron-nickel  alloys  ai  1600  °C.|60Win) 


(2.23) 
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Kp  = a„  - Co  - exp(  - AG,°/RT  ) 
solubility,  C„„,  in  Ihe  dilute  solution  yields 


where  AH,0  and  AS,°  are  the  enthalpy  and  the  entropy  for  oxide  formation  from 
saturated  solid  solution,  respectively.  This  equation  shows  Ihe  temperature 


dependence  of  terminal  solubility. 


refractory  metals  and  their  alloys  containing  oxygen  impurity  are  so  low  that  it 
should  generally  be  expected  that  in  ordinary  vacuum  a pick-up  of  oxygen  from 


residual  gas  atmosphere  occurs  rather  than  a removal.  Therefore,  in  many  cases,  the 
deoxidation  of  refractory  metals  are  performed  via  ihe  reduction  by  carbon  or 
hydrogen,  and  the  suboxide  vaporization. 


The  deoxidation  of  refractory  metals  can  be  achieved  by  the  use 
which,  particularly  under  reduced  pressures,  has  a very  high  affinity  for  ox 


[C]„  + [O],,  = CO(g)  (2.25) 
where  [C]M  and  [0)M  arc,  respectively,  carbon  and  oxygen  dissolved  in  the  metal. 
Fig.  2.6  shows  the  vapor  pressures  of  pure  Nb,  Nb-O,  and  Nb-C-0  sold  solutions 
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Oxygen  in  ppm 


Equilibrium  vapor  pressure  over  (he  Nb-0  and  Nb-C-0  solid 
solutions  at  2273  K.JSOOno]  Vapor  pressure  lines  arc  taken 
after  line(l)|68Fro2J,  line(2)[77Ucd],  Iines<3)  and  (4)[64Gcb|, 
linc(5)|59Spe],  and  !ine(6)[36Rei). 
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wilh  up  10  1 al%  0 ai  2273  K.  A comparison  of  vapor  pressures  of  Nb-0  and  Nb- 
C-O  solid  solutions  indicates  that  CO  pressure  in  equilibrium  with  Nb-C-0  solution 
is  roughly  two  orders  of  magnitude  greater  than  NbO,  vapor  pressure  at  the  same 
oxygen  level.  This  implies  that  deoxidation  due  to  CO  degassing  is  much  more 
effective  than  that  due  to  niobium  oxide  evaporation  which  will  be  described  later. 

On  the  other  hand,  the  deoxidation  by  hydrogen  is  done  according  to  the 
reaction  following  equation. 

[0|„  + H,(g)  = H,0(g)  (2.26) 

with  the  equilibrium  constant,  K, 

K = P,bo/(ao  P,b)  . (2.27) 

This  reduction  process  is  of  particular  interest  for  cases  where  the  reduction  by 
carbon  cannot  be  carried  out  because  of  the  formation  of  stable  carbide.  It  can  be 
seen  that  the  higher  the  ratio  of  P„/Plyo  is,  the  more  effective  the  deoxidation  is. 

In  some  metal-oxygen  systems,  high  vapor  pressure  of  oxides  can  be  utilized 
successfully  for  deoxidadon.  One  of  those  systems  is  Nb-0  system.  In  vacuum 
degassing  of  Nb  containing  oxygen  impurity,  the  low  equilibrium  oxygen  partial 
pressure  of  solution  eliminates  the  possibility  of  oxygen  removal  through  the  direct 
oxygen  molecule  degassing.  Instead,  the  evaporation  of  Nb  suboxides  may  lead  to 
the  decrease  of  oxygen  level.  When  niobium  is  annealed  in  vacuum,  some  of  the 
oxygen  striking  the  surface  of  Nb  evaporates  back  into  vacuum  in  the  form  of 
suboxides  while  some  others  are  dissolved  in  Nb.  Since  the  oxygen  content  in 
niobium  increases  with  time  at  a given  partial  pressure  of  oxygen  and  desorption 
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reached  when  Ihc  nnc  of  oxygen  removal  due  to  oxide  evaporation  is  equal  to  the 
rate  of  oxygen  dissolution.  At  this  condition,  the  oxygen  concentration  remains 
constant  at  a level  corresponding  to  the  given  oxygen  partial  pressure  and 
temperature.  Oxygen  concentration  in  solid  Nb  at  steady  state  is  reported  to  be 
lower  as  temperature  increases  or  oxygen  partial  pressure  dccreasesI68Horl|. 

The  possibility  of  deoxidation  via  volatile  oxide  evaporation  can  usually  be 
examined  using  the  evaporation  ratio  suggested  by  Brewer  and  Rosenblatt(62Bre), 
who  developed  the  idea  based  upon  the  thcrmochemical  stability  of  oxide.  The 
definition  of  the  ratio  is  given  below. 


R _ (I  y PMxOy  /X  x Pm»o,L  . 

< Nr/fU— 


(2.28) 


where  PMxOy  is  the  equilibrium  partial  pressure  of  volatile  oxides  being  considered 


and  N„  and  N„  the  molar  traction  of  oxygen  and  metal  in  the  condensed  phase.  For 
favorable  purification,  it  is  recommended  that  R should  be  at  least  10  or  greater. 
Table  3 [820no]  summarizes  the  equilibrium  thermodynamic  data  concerning 
evaporatron  dcoxtdation  and  the  calculated  evaporation  ratio,  R,  for  several  elements. 
The  elements  in  group  IV  have  a large  solubility  and  lose  oxygen  as  an  MO  species 
so  that  the  value  of  R is  small.  The  elements  in  group  V lose  oxygen  as  suboxides. 
Purification  via  suboxide  vaporization  could  be  practical  for  these  elements  except 
V-O  system,  where  the  vapor  pressure  of  the  base  metal  is  relatively  large. 
Chromium  in  group  VI  is  not  a good  candidate  because  of  the  high  vapor  pressure 
of  Cr  itself.  The  refractory  metals  Mo  and  W of  the  group  Vi  have  extremely  high 
R values,  indicating  vety  effective  purification  by  vaporization.  The  major 
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vaporization  oxide  species  are  MO,  type  molecules  and  its  dimer.  These  metals  arc 
characterized  by  the  small  solubility  and  low  activity  of  oxygen. 


The  Nb-N  phase  diagram  is  shown  in  Fig.  2.7.  Up  to  60  at%  N,  there  are 
two  nitride  phases,  Nb,N,  and  NbN.  Nb,N  has  a hexagonal  close  packed  structure 
and  homogeneity  range  of  28.5  to  32.4  at%  N.  The  NbN  phase  undergoes  a 
transformation  from  face  centered  cubic  to  hexagonal  close  packed  at  approximately 
1300  °C.  The  solubility  of  nitrogen  in  Nb,  C„.  is  given  as  a function  of  partial 
pressure  of  nitrogen  and  temperature  by  the  following  equations.  For  nitrogen 
solubility  in  a -Nb[68Ftol),  it  is 

log  Cy  = 1/2  log  Pw  -3.10  + 9300/T  <1770<T<2470  K)  (2.29a) 
and  for  N in  liquid  Nb|73Rev],  it  is  given  as 

log  Q,  = 1/2  log  P„,  + 0.28  (<  15  at%,  T = 2800  K)  (2.29b) 
with  Py,  in  mbar,  T in  K and  C„  in  at%.  Some  thermodynamic  data  associated  with 
Nb-N  system  is  summarized  in  Table  4. 

The  pressure-concentration  isotherms  are  shown  in  Fig.  2.8  in  which  a slope 
of  2 for  the  region  of  a-solid  solution  implies  that  nitrogen  is  dissolved  in  atomic 
form  in  the  solid  solution.  The  kink  of  the  isotherms  in  the  lower  regime  of  G, 
gives  the  concentration  at  which  niobium  is  saturated  with  nitrogen  at  the  given 
temperature  and  also  the  subnitridc  phase,  Nb,N,  starts  to  form.  For  constant  N, 


pressures. 
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Fig.  2.7  Nb 


- N phase  diagram.I81Schl) 
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Table  4.  Themiodynamic  data  associated  with  Nb  - N systen 


Equilibria  AG° 

C in  at%.  P in  loir,  T in  K in  cal 

References 

1/21N.)  - INI. 

logC*  = l/21ogP,  - 3.1  + 930U/T  -42500  + 16.8T 

( 1730  - 2470  K ) 

[68FroIl 

*logC„  = l/21ogP»  + 0.22 
( T = 2800  K;  < 15  at%  ) 

2gNb>  + l/2(N,l  ° <Nb,N>  (1730  - 2470  K 1 

(73Rev  | 

logPN,  = 12.6  - 2840071'  -65100  + 22.2T 

[68Fcoll 

log  C*  = 3.17  - 4920/T  +22500  - 5.35T 

[68Froll 

* Solubility  in  liquid  phase  and  P in  mbar. 

Fig.  2.8  Pressure  - 


- N sysiem.{66Geb) 
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temperature.  The  other  kink  in  each  isotherm  indicates  the  concentration  at  which 
nitrogen  starts  to  dissolve  in  Nb,N  in  the  form  of  solid  solution.  In  Fig.  2.9.  the 
solubilities  in  niobium  at  2050  ”C,  just  below  and  above  its  melting  point,  arc 
plotted  against  the  N,  partial  pressure.  It  is  seen  that  at  the  given  N,  pressure  more 
amount  of  nitrogen  is  dissolved  in  the  melt  than  in  solid. 

Nb  - O System 

The  Nb  - O phase  diagram  is  shown  in  Fig.  2.10.  Three  different  oxides 
exist.  NbO,  NbO,,  and  N'b.Ov  The  melting  points  of  the  oxides  are  1945.  1925.  and 
1495  °C,  respectively.  The  maximum  solubility  of  oxygen  in  a-Nb  is  7.5  at.%  at 
1915  °C  The  O,  equilibrium  pressure  over  the  a solid  solution  is  so  low  that  direct 
measurements  of  the  solubility  data  cannot  be  performed  even  under  high  vacuum. 
Instead,  the  equilibrium  data  can  be  obtained  from  the  temperature  dependence  of 
the  terminal  O solubility  and  AG°  values  of  oxide  formation.  The  solubility  in  a- 
solid  Nb  is  given  by  equation[68Fro!]. 

log  Co  = 1/2  log  P„  - 4.45  + 20020/T  (2.30) 

( 1570  < T < 2270  K ).  with  P„,  in  mbar.  T in  K and  C„  in  at.  %. 

The  terminal  solubility  C,.„  in  a-Nb  is  given  as)68Frol| 

C._  = 1.67  - 1680  T,  (1070  < T < 1770  K)  (2.31) 

The  vapor  pressure  of  the  volatile  oxides  over  Nb  containing  O were  determined  to 
bc|68Frol] 

NbO(g):  log  PSV)  = log  Co  + 9.12  - 2840OT  (2.32a) 


(1870  < T<  2470  K), 
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Fig.  2.9  Nb  - N systi 
Nb.[81SchI] 


n.  log  I'  - log  Cs  diagram 


melting  point  of 
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Fig.  2.10  Nb  - O phase  diagram.[60Ell,  68Frol] 
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NbO,(g)  : log  P**,  = 2 log  Q,  . „ + 7.12  - 24300/T  (2.32b) 

(1870  < T<  2470  K) 

Thermodynamic  data  concerning  Nb-O  system  are  summarized  in  Table  5.  The 
pressure  - concentration  isotherms  for  the  steady  states  of  niobium  annealed  in  03  is 
plotted  logarithmically  in  Fig.  2.11|77Sch].  The  isotherm  shows  a slope  of  1 at 
constant  temperature,  Co  « Po;.  Evaporating  species  consists  of  two  oxides,  NbO 
and  NbOj.  and  NbO  is  a predominant  at  lower  concentration  and  higher 
temperatures.  The  percentage  of  NbO  is  indicated  by  dot-dashed  lines.  The  existence 


of  the  a-solid  solution  limited  by  the 
in  the  upper  part  of  the  graphs). 


of  steady  state  is  restricta 
solvus  and  solidus  lines  ( t 
Furthermore,  at  low  temperatures  the  steady  state  conditions  cannot  be  realized  any 
longer  due  to  the  negligible  oxide  evaporation,  leading  to  a continuous  oxygen 
absotption  and  oxide  scale  formation.  Only  at  very  low  O2  pressure  is  a steady  state 
still  available,  but  due  to  the  extremely  low  reaction  rate,  very  long  period  is 
necessary.  The  oxygen  concentration  in  the  steady  state  condition  is  found  to  be 
related  to  the  oxygen  partial  pressure  and  temperature  as  follows(77Schl. 


Co  = Pa  • 9.1  • 10"  exp(-  502000/RT  ) 

( 2170  < T < 2540  K ) 

where  C„  is  in  al.%  and  P„  in  mbar,  and  R is  8.31  J/mol.K. 

treatments.  On  the  other  hand,  it  should  be  mentioned  that,  in 
equilibrium,  an  additional  metal  loss  always  appears  in  steady  si 


(2.33) 


constant  oxide  evaporation  rate,  which  can  be  larger  by  several  orders  of  magnitude 


Tabic  5.  Thcimodyn 


Equilibria  AG°  References 

C in  al%,  P in  lore.  T in  K in  cal 


1/2(0.)  = IOI. 

logCo  = 1/2  logP,,  - 4.5  + 20200/T  -92300  + 23.0T  [68Frol] 

1/2(0)  = IOI„ 

logCo  = 1/2  logP„  - 9.4  + 33400/T  -152900  + 39.0T  [68F10II 

(WbO)  = 101.  + <;NbP 

logQ,  = logPsK>  - 9.0  + 28400/T  -13000  + 37.1T  [64Frol] 

l/2(NbQ,)  ° 101..  + l/2<Wb> 

logCo  - l/21ogP™.  - 3.5  + 12150/T  -55500  + 18.5T  [68F10I] 

<Nb>  + 1/2(0.)  = <NbO> 

logP„  = 12.3  - 43700/T  -100000  + 21.5T  [68F10I] 

<Nb>  + (O)  = <NbO> 

logP0  = 11.1  - 35100/T  -160600  + 37.5T  [68F10II 

<NbO>  = <Nb>  + 1Q1„ 

logCo  = 1.67  - 1680/T  +7700  + 1.5T  I72FrolJ 

Sreadv  stale 

O,  logCo  = logP„,  - 3.35  + 16700/T  [69Frol] 

( ) gas  phase,  [ ] imersiiiial  in  solid  solution,  and  < > solid  phase. 


Oxygen  pressure  in  Torr 


Fig.  2.11  Steady  state  during  annealing  of  Nb  in  O,  atraospherc.|77Sch| 
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than  metal  evaporation  rale,  as  can  be  seen  from  Fig.  2.12[69Fro,  70Frol],  Tile 
evaporation  rale  is  given  by  the  thickness  decrease  of  base  metal  per  hour  in  }im/h, 
It  can  be  seen  that,  in  the  pressure-temperature  range  of  steady  state,  NbO  and 
NbOj  gases  arc  predominant  evaporating  species.  However,  at  very  high  temperature 
and  the  low  pressure  region,  it  is  possible  for  the  pure  metal  evaporation  to  exceed 
the  oxide  evaporation. 

Nh  - C system 

The  Nb-  C phase  diagram  is  shown  in  Fig.  2.13.  There  exists  two  carbide 
phases,  namely,  hexagonal  close  packed  Nb,C  and  NaCI-type  cubic  NbN.  The 
maximum  solubility  of  carbon  in  a-Nb  is  3.5  at.%  at  about  2300  °C.  In  the  Nb-C 
system  only  the  terminal  solubility  of  the  a solid  solution  can  be  determined  as  a 
function  of  temperature  as  follows: 

log  Cc  = 3.65  - 7600/T  (1770  < T < 2470  K)  . (2.34) 

During  annealing  of  niobium  in  CO  atmospheres,  thermodynamic  equilibrium  is 
established  according  to  the  reaction. 

(CO),  = [QU+  [C]„.  (2.35) 

The  CO  pressure  over  Nb-C-O  solid  solutions  are  given  by  the  following 
expression. |6.SFro2) 


log  Pro  = log  Cc  + log  C„  + 6.32  - 14700/T  (2.36) 

(1870  < T < 2170  K),  with  Pc  in  mbar  and  Q-  and  Q,  in  at%. 

Other  thermodynamic  data  on  Nb-C  system  is  given  in  Table  6. 

In  Fig.  2.14,  experimentally  determined  CO  equilibrium  pressure  isotherms 
for  Nb-C-O  solid  solution  with  the  concentration  ratio  C/O  =1  arc  graphically  shown 
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of  Nb  in  0,.[69Fro) 


18  22 
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Fig.  2.13  Nb  - C phase  diagram.|81Sch] 


able  6.  Thermodynamic  data 


Equilibria  AG°  References 

C in  at%,  P in  ioit,  T in  K in  cal 


SG>  - I.CI„ 

2<Nb>  + <C>  = <;NbiC> 

<Nb.C>  = 2<Nh>  + in, 
logC,  = 3.65  - 7600/T 


1/210,1  + <C>  = (CO) 


-3100  - 9.5T 


+34800  - 7.53T 
• 14700/T  -67000  + 33.5T 
-28400  - 20.0T 


I64Gebl 
[68Froll 
[66Geb] 
[68  Fro  1 1 
[68Froi] 


phase. 


Fig.  2.14  Equilibrium  pressure  of  CO  over  Nb-C-O  solid  solution  with  the 
of  C to  O equal  to  unity.[76Fro) 


logarithmic  scale.  Equilibrium  C or  O 


of  CO  decreases  with  increasing  temperature,  indicating  that  CO  dissolution  in  Nb  is 
an  exothermic  reaction.  The  CO  equilibrium  pressure  over  solid  solution  of  C and  O 
containing  a fraction  of  at%  at  2000  to  2600  K is  relatively  high,  malting  it  possible 
for  reaction(2.35)  to  be  used  for  removing  carbon  from  niobium  by  a high  vacuum 
treatment.  For  a given  partial  pressure  of  CO,  C.-C.  is  constant  at  a constant 
temperature.  Therefore,  carbon  can  be  removed  by  annealing  the  specimen  in  O, 
atmosphere. 


The  rate  of  reaction  is  formally  d 
component  per  unit  volume  with  respect  to 


id  as  the  change  in  moles  of  a 
e.  The  mathematical  expression  for 


r=  ~v  ~ar  - ~ar 

where  V.  N.  and  C are  the  volume,  the  num 
concentration  of  the  species  involved  in  the  re 
coefficients  for  two  reactants  are  different,  the  n 


(2.37) 

t of  moles,  und  the  molar 
ition.  If  the  stoichiometric 
i expressed  in  terms  of  one 
i terms  of  the  other.  Consider 


changes  for  the  four  : 


) (2.38) 

and  products  are  related  as  follows. 


1 dID] 


(2.39) 
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where  [A]  is  the  conceniration  of  reactant  A.  etc. 

In  Older  for  reaction  to  take  place,  each  reactant  should  be  close  enough  to 
make  collision.  Therefore,  the  reaction  will  be  proportional  to  the  frequency  of 
collision  which,  in  turn,  is  associated  with  the  concentration  of  each  species  in  the 
following  way: 


I d[A] 

TTT 


k [Af  [B]- 


(2.40) 


where  n and  m are  the  order  of  reaction  with  respect  to  |A]  and  [B|  , respectively, 
and  k is  the  rate  constant.  The  sum  of  n and  m is  called  the  order  of  the  overall 
reaction.  The  order  and  the  stoichiometric  coefficients  in  the  reaction  equation  are 
not  necessarily  intenelated,  i.e.,  it  is  not  required  that  n=a  and  m=b.  In  general,  the 
order  is  determined  experimentally.  Since  concentration  data  is,  in  practice, 
measured  as  a function  of  time,  the  integral  form,  rather  than  the  differential  form, 
of  reaction  rate  equation  is  more  frequently  used  in  determining  the  order  of 
reaction  because  data  can  be  plugged  into  it  directly. 


Rate  Controlling  Step 

As  shown  in  many  kinetics  studies,  the  formation  of  the  final  products  from 
the  original  reactants  usually  occurs  in  a series  of  relatively  simple  steps  which  may 
be  physical  or  chemical  in  character.  In  order  to  understand  the  gas-metal  reaction 
in  terms  of  kinetics,  it  is  important  to  consider  the  reaction  mechanism  which  is  the 
sequence  of  steps  that  describe  how  the  final  products  are  formed  from  the  original 
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reactants.  The  rale  of  the  individual  steps  will  normally  be  different,  and  [he  rate  of 
the  overall  reaction  will  be  determined  primarily  by  the  slowest  among  these  steps. 
Herein  lies  the  explanation  for  the  difference  between  the  order  of  reaction  and  the 
stoichiometric  coefficients.  Together  with  the  reaction  mechanism,  it  is  also 
important  to  understand  the  concept  of  rate  controlling  step.  While  the  notion  of 
"bottle  neck"  process  in  a series  of  subprocess  is  intuitively  obvious,  it  is  often 
carelessly  applied  or  stated.  Since  a clear  understanding  of  the  implication  of  a rate 
controlling  or  the  slowest  step  is  central  to  the  correct  interpretation  of  coupled 
physical  and  chemical  reaction  phenomena,  the  concept  of  rate  controlling  step 
motivated  by  C Wagner|70Wagl  will  be  mentioned  below. 

Consider  a sequence  of  consecutive  steps  of  a reaction  involving  species,  A 
as  die  reactant,  Z„  Zn*  a * as  the  intermediates,  and  species  B as  the  product. 

A = Z,  = Zj  = ««*  = B (2.41) 

where  steps  are  numbered  by  1,  II,  • • • ,n.  Since  the  concentrations  of  the 

intermediates  are  small  in  comparison  to  those  to  A and  B,  after  a relatively  short 

induction  period,  a steady  stale  with  virtually  equal  net  rate,  v„  v„  • • • of  the 

individual  steps  is  established.  Thus,  the  net  rate,  v,  of  the  overall  reaction  A = B 

V “ V,  = v„ = v.  (2.42) 

The  above  statement  is  not  consistent  with  the  general  definition  of  rate  determining 
step  and  therefore  must  be  modified.  This  is  possible  by  introducing  the  concept  of 
the  virtual  maximum  rate,  v,'  , of  each  individual  step.  The  virtual  maximum  rate  is 
defined  as  the  rate  which  would  be  observed  if  for  all  steps  other  than  step  i virtual 


equilibrium  were  established,  e.g.  with  the  help  of  appropriate  catalysts  (and  all 
other  irreversible  steps  were  comparatively  fast).  Then  the  concentration  Q,'  of 
species  Z,.  in  the  case  of  equilibrium  for  ail  steps  preceding  step  i equals  to 

Cu*  - C.  (K,  K„  • • • K,,)  (2.43) 

where  K,  denotes  the  equilibrium  constant  of  step  i.  Furthermore,  the  concentration 
of  species  Z in  the  case  of  equilibrium  for  all  steps  following  step  i is  given  by 

C;  = C.  (K,„K.,  • • • KJ-'  (2.44) 

The  equilibrium  constant  K of  the  overall  reaction,  A = B,  is  related  to  that  of  each 
step  as  follows: 

K = (K,  K„  • • • KJ  (2.45) 

For  each  individual  step,  virtual  maximum  rate,  v,\  is  then  given  by 
v,'  = k,CA  - k,'C,‘  = k,  ( C4  - C,7K, ) = k,(C,  - CVK) 
v„‘  = k„C,"  - kn'C,'  = k„  ( C,-  - C,'  /K„)  = k„K,(CA  - CJK)  (2.46) 


v.'  = k.0,,-  - k.’C  = k,  ( C,‘  - C.-/K.  ) = k,KfclK.a*"K,(CA  - CVK) 
where  k,  and  k,'(=  VK)  indicate  the  rate  constants  of  the  forwaid  and  the  backward 
reaction  of  step  i.  respectively.  On  the  other  hand,  if  equilibrium  for  the  individual 
steps  is  not  established,  it  yields 

v = k,  ( CA  - CVK,  ) 
v = k„  ( C,  - CVK„  ) 


(2.47) 


k,  ( C.,  - CJK.) 


Dividing  through  corresponding  side  of  equation  (2.46)  and  (2.47)  and  adding 
die  individual  quotients,  it  yields 


(2.48) 


v " ' i/v/  + i/v‘- . . . + i/v; 

This  equation  implies  that  the  observed  reaction  rate, 
the  individual  rates,  v,\  and  is  practically  equal  to  the 
step  j if  v,'  is  much  less  than  the  virtual  maximum  rates 


(2.49) 

v,  inherits  characteristics  of 
of  all  other  steps. 


v = v/  if  »;  « v„'  (2.50) 

Therefore,  it  can  be  justified  that  the  overall  reaction  rate,  v,  depends  only  on  the 
parameters  governing  vj.  becoming  insensitive  to  the  characteristics  of  all  of  the 
other  rate  equations.  v,'(i*j). 


In  deriving  die  rate  equation,  dC/dt  = k[Aj,  the  effect  of  concentration 
difference  between  the  initial  and  the  final  states  is  directly  incorporated  into  the 
concentration  terms  as  a driving  force  for  the  reaction,  but  the  influence  of  the  other 
parameters  is  represented  in  the  rate  constant,  k.  Such  parameters  include  the 
resistance  to  mass  transfer  in  the  gas  phase  or  liquid  phase,  temperature,  physical 
conditions  of  specimen  such  as  surface  to  volume  ratio  and  so  on.  The  effect  of 
barrier  to  the  mass  transport  will  be  discussed  later  in  detail.  Temperature  is  one  of 


the  most  important  variables  to  be  included  in  the  rate  constant,  particularly,  for  the 
intcrfacial  elementary  steps. 

Temperature  dependency  of  rate  constant,  k,  for  an  elementary  process  obeys 
the  Arrhenius  equation.  Thus,  it  is  given  by 

k = k„  exp(-Q/RT)  (2.51) 

where  ko  is  a prccxponential  constant. 

Q is  the  activation  energy. 

R is  the  gas  constant. 

T is  the  absolute  temperature. 

Combining  cquation(2.51)  and  (2.40)  yields 

-dIAl/dt  = k,  exp(- QflO)  [Af[Br  (2.52) 

it  should  be  noted  that  this  expression  is  rigorously  limited  to  an  elementary  process 
because  the  Arrhenius  equation  is  so  restricted.  However,  the  exponential  effect  of 
temperature  often  accurately  represents  experimental  rate  data  for  an  overall  reacdon, 
even  though  the  activation  energy  is  not  clearly  defined  and  may  be  a combination 
of  Q values  for  several  elementary  steps. 

The  different  value  of  aedvadon  energy  for  each  elementary  step  leads  to  a 
different  temperature  dependence  of  each  step,  resulting  in  the  switch  of  rate 
controlling  step  from  one  to  another  in  different  temperature  regime.  According  to 
Fig.  2.15,  the  elemcntaty  step  B with  higher  activation  energy  may  be  the  slowest  at 
low  temperaturesfshown  by  range  a in  the  figure)  and  could  well  become  rate 
determining  for  the  total  degassing  reaction.  On  the  other  hand,  at  high 
tcmpcralures(range  c)  the  step  C with  lower  activation  energy  is  the  rate 


Fig.  2.15  Rale  of  a rcaciion(A)  consisting  of  two  partial 
C)  which  are  both  rate  delermining|71Winl]. 


determining.  In  the  intermediate  temperature  range,  b,  both  elementary  steps  proceed 
at  about  the  same  rate  so  that  the  whole  reaction  process  becomes  under  control  of 
both  steps  in  a mixed  mode. 


In  the  earlier  section,  it  was  discussed  how  the  overall  reaction  is  related  to 
each  elementary  step  and  the  overall  reaction  rate  is  associated  with  the  individual 
substep  reaction  rates.  In  the  following  section,  the  relationship  between  the  overall 
reaction  and  individual  elementary  steps  will  be  discussed  from  the  energetic  point 
of  view.  For  example,  the  energy  diagrams  for  A,  dissolution  reaction  are 
schematically  illustrated  in  Fig.  2.16  where  the  energy  of  one  mole  A,  or  two  g- 
atoms  A is  plotted  against  the  distance  of  the  molecules  or  atoms  from  the  metal 
surface  for  two  types  of  metal-  gas  systems.  Fig.  2.16(a)  shows  the  energy  change 
for  a system  with  a large  negative  heat  of  solution  AH°  such  as  Nb-N,  Ta-N,  and 
V-N  (Vb  group  - N),  whereas  Fig,  2.16(b)  is  for  a system  with  a large  positive 
value  of  AH0  such  as  Cr-N,  Mo-N,  and  W-N  (Vlb  group  - N).  During  gas 
dissolution,  an  A;  molecule  passes  at  first  through  a shallow  potential  minimum  of 
molecular  adsorption  and  reach  the  activated  state  for  dissolution  as  A atoms.  The  A 
atoms  are  then  chemisorbed,  penetrate  the  metal  surface,  and  diffuse  in  the  bulk  of 
the  metal.  During  desorption  the  same  steps  are  passed  in  the  reverse  direction.  Q„ 

and  Q„.  denote  the  apparent  activation  energies  for  absorption  and  desorption.  O 

Q,',  Q„  and  Q*,  arc  the  true  activation  energies  for  dissociative  chemisorption. 


Fig-  2-16  Energy  diagrams  for  reversible  absorption  and  desorption  of 
diatomic  gases  A,.  AH„  dissociation  energy,  AH_„  AH,  and 
AH1C|  heat  of  chcmisotption,  solution,  and  segregation;  Q*  and 
O.  experimental  activation  energy  of  absorption  and  desorption; 
0- — Q„  or  0,.  Q, true  activation  energy  of  dissociative 
chemisorption,  transfer  through  metal  surface,  and  diffusion. 
|80Fro| 


(a)  exothermic  solution  reaction,  (b)  endothermic  solution 
reaction. 
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transfer  through  the  metal  surface  in  absorption  and  desorption,  and  diffusion  of  A 
atom,  respectively.  AH”  indicates  the  heat  of  solution. 

The  differences  in  reaction  rates  and  activation  energies  for  nitrogen 
absorption  and  desorption  between  group  Vb  metal  and  group  VIb  metals  can  be 
correlated  with  the  different  sign  for  the  heat  of  solution  AH°  of  each  group.  In  the 
case  of  Nb,  AH0  has  a large  negative  value,  i.e.,  the  solution  of  nitrogen  is 
exothermic.  So,  if  the  shallow  potential  minimum  of  weak  adsorption  is  neglected, 
the  apparent  activation  energies  for  nitrogen  absorption(Q„)  and  dcsorptionfQ*.)  are 
related  to  that  for  elementary  steps  and  the  heat  of  soIution(AH°)  as  follows: 

Q.  = Q_  (2.53) 

Q*.  “ 1 AH0  1 + Q*  . (2,54) 

With  these  energetic  relationships,  it  can  be  predicted  easily  that  desorption  rate  is 
slower  than  absorption  rate  in  the  system  with  exothermic  solution  reaction  if  the 
driving  force,  i.e,  concentration  differences  between  the  initial  and  the  final  state,  is 
identical.  In  contrast,  when  AH0  is  positive  such  as  for  Mo  and  W,  i.e.,  the  solution 
reaction  is  endothermic,  those  are 

Q»  = 1/2  AH°  + Q„  (2.55) 


with  Q,  assumed  to  be  equal  to  In  the  previous  equations,  AH0  refers  to  the 
value  for  one  mole  of  A,  and  thus  numerical  values  are  twice  the  values  given  in 


SI 

unil  of  cal/g-atoms.  Based  upon  (he  same  reasoning  as  for  exothermic  reaction, 


aspects  of  high  temperature  gas/metal  re 


d rate  laws.  Most 


rate.  These  are: 


Hence,  step(i)  and  step(iii)  are  physical  in  character  while  step(ii)  is  chemical.  The 
composition  profiles  in  the  gas  phase  and  metal  are  shown  schematically  in  Fig. 


ather  refractory  metals. 
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Fig.  2.17  Concentration  and  partial  pressure 
in  degassing  process. 


profile  for  different  types  of  : 
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Mass-Transfer  in  the  Gas  Phase 

When  species  A in  a gas  stream  of  A + B mixture  is  under  the  influence  of 
concentration  gradient,  the  molar  flux  of  A relative  to  stationary  coordinate  is  made 
up  of  two  ports:  -CDa,VXa  which  is  the  diffusion  contribution  and  XA(  N,  + N„) 
which  is  the  molar  flux  resulting  from  the  bulk  motion 


Na  = - CD„VXa  + XA(  N.  + N„  ) (2.57) 

where  Na  and  NB  are  molar  fluxes  of  A and  B,  C is  the  concentration  of  A + B. 

is  the  dilfusivity  of  A in  A + B mixture,  and  XA  is  the  molar  fraction  of  A.  If 
bulk  motion  resulting  from  volume  change  caused  by  surface  chemical  reaction  docs 


not  exist,  the  second 
of  diatomic  gases  like 
gas  boundary  layer  car 


term  in  equation  (2.57)  ca 
N,  or  H,  belongs  to  this 
r be  described  as  follows. 


in  be  neglected.  Dissolution  reaction 
category.  Then  mass  transfer  in  the 


Na  = -CDa„VXa  (2.58) 

Assuming  that  the  A + B gas  mixture  behaves  ideally,  integration  of  equation  (2.58) 
over  the  boundary  layer,  8,  gives 


where  b and  i refer  to  the  bulk  and  interface,  respectively,  and  7',  is  equal  to  l/2(Tb 
+ T).  This  equation  can  be  rewritten  in  terms  of  wt%/dt  and  partial  pressures  as 


where  CA  is  the  concentration  of  A in  the  gas  in  wt%,  MA  is  the  atomic  weight  of 
A,  A is  the  surface  area,  W is  the  weight  of  the  sample,  and  PA  is  the  partial 
pressure  of  A.  Since  CA’  KPW  equation  (2.60)  yields 


dC^di  = 


(2.61) 


This  equation  implies  that  if  gas  phase  transport  is  rate  controlling  in  diatomic  gas 


dissolution  reaction,  then  the  reaction  should  be  second  order  with  respect  to  CA. 

On  the  other  hand,  the  molecular  flux  due  to  bulk  motion  should  be  taken 
into  account  if  the  surface  reaction  results  in  volume  change.  For  example,  consider 

the  decarburization  by  oxygen;  C + O,  = 2 CO.  In  this  reaction,  the  net  flux  of 

gas  takes  place  away  from  the  metal  since  two  moles  of  CO  gas  are  generated  for 
each  O,  molecule  reacting  with  the  metal.  Thus,  the  net  flux  away  from  the  metal 
should  decrease  the  rate  of  transport  of  reacting  species.  In  such  case,  equation 
(2.57)  can  be  applied  without  neglecting  the  second  term. 


where  C is  the  concentration  of  O,  + CO,  Do,«,  is  die  diffusivity  of  O,  in  O,  + CO 
mixture,  and  is  the  molar  fraction  of  O,,  N„  and  Nco  are  molar  fluxes  of  O, 
and  CO,  respectively.  Plugging  Nn,  = - 2N„,  in  integrating  the  equation  (2.62)  over 
the  gas  boundary  layer,  6,  it  gives 


where  P„  is  the  partial  pressure  of  oxygen.  Using  ln(l+P„/P)  - P„,  when  Po,  is 
small  compared  with  total  pressure  of  1 amt,  equation  (2.63)  can  be  reduced  to 


No.  = -CDo^odXo/dx  + X„<  N„  + Nco) 


(2.62) 


RT, 


-D°»"’  (Po,*  -Po,') 


8 


(2.64a) 


(2.64b) 
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-^V  4c  <«-«> 

On  the  other  hand,  close  examination  of  equations(2.61)  or  (2.64b)  indicates 
that  only  D„  and  5 ate  unknown  parameters.  An  estimation  of  gas  phase  mass 
transfer  coefficient,  k„,  which  is  defined  as  D„/8,  can  be  made  from  the  Steinbcrger 
and  Treybal  equation  (60Stel,  which  relates  the  Sherwood  number  to  other 


Sh  = 2 + 0.5  (Gr'Sc)“  + 0.35(RcSc“!y'“ 
and  Gr'  in  the  above  equation  is  defined  as 


Gr-  = GrM  + (Sc/Pr)05  GrH 


GrM;  Grasf 
GrH;  Grash 
Pr,  Prandtl 


her,  (=  M/D) 

iber  for  mass  transfer  (=gd’p<(p,-pJ/P’) 
ber  for  heat  transfer  (=gd’p,(TrT>)/pf,T,) 

(=pcyk), 


Re;  Sphere  Reynolds  number  (=p^V/p,) 

Sc;  Schmidt  number  (=p/pD), 

M gas  phase  mass  transfer  coefficient  (cm  s ') 


d;  diameter  of  specimen  (cm) 


(2.65) 

(2.66) 


cm1  s’1) 
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ji;  viscosity  of  gas  (g  cm-1  s'1) 

T;  temperature  (K) 

V;  velocity  of  gas  (cm  s ') 

C,;  heat  capacity  of  gas  (J  g'1  K ‘) 
k;  themial  conductivity  of  gas  (J  cur’  s'1  K1) 


g;  acceleration  due  to  gravity  (cm  s‘) 
i;  interface 
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tested  by  a number  of  high  temperature  expcrimcnts[67Bak,  71Fru,  85Lce,  85Rao, 
86Bat|.  For  example,  in  the  deoxidation  study  of  copper  by  !%  hydrogen  in 
nitrogcn[71Ric]  it  has  been  found  that  at  temperatures  between  1400  and  1600  “C 
the  values  of  Iq,  derived  from  the  experiments  agree  well  with  the  values  derived 
from  the  equation(2.65).  Good  agreement  between  the  experimental  and  calculated 
values  for  various  flow  rates  is  also  repotted  in  the  copper  oxidation  at  1600  °C 
under  flic  nitrogen  atmosphere  containing  a small  amount  of  oxygen.|70Gle  ,71Ric] 

Therefore,  in  the  light  of  all  the  above  discussion  it  can  be  summarized  that 
if  gas  phase  transport  is  the  rate  controlling  step,  then 

(1)  absorption  (or  desorption)  rate  would  be  proportional  to  (Pf  - PA),  or 
(CJ-C*)  in  the  dissolution  reaction  of  diatomic  gas,  i,c„  second  order 
with  respect  to  [CJ. 

(2)  gas/metal  reactions  can  be  assumed  to  be  controlled  by  transport  in  the 

gas  phase  when  the  rate  is  markedly  influenced  by  gas  velocity  because 
ko  is  related  to  Reynolds  number  which  includes  the  gas  velocity  term. 

The  mass  transport  in  the  interior  of  the  liquid  phase  is  affected  by 
convection,  which  may  be  generated  by  artificial  means  such  as  induction  stirring  or 
blowing  gases  through  the  melt.  It  may  also  be  induced  spontaneously  by 
differences  in  the  density  of  the  melt  between  the  surface  and  the  interior,  which  is 
caused  by  temperature  difference  resulting  from  heat  loss  from  the  surface  of  the 
melt.  The  convection  can  accelerate  equalization  of  tile  concentration  within  the 
melt.  But  it  is  to  be  realized  that  the  most  important  process  in  liquid  phase 


transport  is  the  transport  to  the  gas/liquid  phase  boundary,  taking  place  by  diffusion 
in  the  liquid  boundary  layer. 

The  boundary  layer  mass  transport  rate  is  proportional  to  the  concentration 
difference  between  the  interior  and  the  surface  of  the  melt.  The  rate  is  therefore 

dCJdt  = - ftA/VXQ.  - CJ  (2.67) 

where  C,  is  the  concentration  in  the  interior  of  the  melt  and  C„  is  the  concentration 
at  the  surface  which  is  in  equilibrium  with  the  partial  pressure  of  the  substance  to 
be  removed  outside  the  melt,  and  p is  the  mass  transfer  coefficient  in  liquid  which 
has  the  dimension  of  Cm/scc.  In  a similar  manner  to  that  in  gas  phase  transport,  p 
is  related  to  diffusion  in  the  boundary  layer  as 

P = D/8  (2.68) 

where  D is  the  diffusivity  of  element  of  concern  in  the  melt  and  8 is  the  thickness 
of  boundary  layer.  Thus  it  can  be  seen  that  the  mass  transfer  coefficient  depends 
mainly  upon  diffusivity  of  a component  of  interest  in  the  liquid  and  hydrodynamic 
factors  affecting  the  thickness  of  boundary  layer.  In  contrast  to  gas  phase  transport, 
there  is  no  reliable  means  to  estimate  the  mass  transfer  coefficient  for  the  interior  of 
the  liquid  metal  because  of  the  poor  knowledge  of  liquid  flow  during 
electromagnetic  levitation.  However,  it  is  well  known  that  a great  deal  of  movement 
goes  on  in  both  the  surface  and  the  interior  of  the  levitated  droplet  although  this 
movement  varies  to  a great  extent  with  the  metal  and  with  the  size  of  liquid 
droplet.  Two  gram  drops  of  iron,  for  example,  have  been  observed  to  pulsate  at  the 
frequency  of  about  40  Hz  with  the  diameter  changing  each  cycle  by  a factor  of 
two.(71Ric] 
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Some  experimental  results  make  it  certain  that  such  movement  enhances  the 
mass  transfer  in  the  liquid  metal.  In  the  experiment  on  deoxidation  of  the  levitated 
copper  droplet  in  hydrogen  and  oxidation  of  nickel  in  oxygcnl7IRicl,  the  mass 
transfer  coefficient  of  0.2  Cm/sec  was  obtained.  This  value  permits  an  8 mm 
diameter  drop  weighing  about  2 gram  to  come  halfway  to  equilibrium  concentration 
in  0.46  second.  In  comparison,  a rigid  drop  of  similar  size  would  bring  its  average 
concentration  halfway  to  equilibrium  in  48  sec  for  a diffusion  coefficient  of  Iff* 
CraVsec.  The  diffusion  coefficients  of  oxygen  in  copper  at  1100  °C  and  in  iron  at 
1550  °C  arc  4- 10’  and  8-10’  Cm’/sec,  respectively.  Thus,  it  can  be  seen  that  fluid 
motion  in  the  levitated  drops  enhances  material  transport  to  surroundings  100  times 
faster  than  the  droplet  without  motion.  Another  interesting  results  of  the  mass 
transfer  in  levitated  melt  is  shown  in  Lee  and  Rao's  experimentsl82Lee]  in  which 
the  effective  diffusivity  of  carbon  in  an  iron-carbon  droplet  levitated  in  the 
electromagnetic  field  was  about  10  times  as  large  as  that  in  the  similar  melt 
contained  in  the  crucible.  This  also  indicates  that  a vigorous  motion  present  in  the 
levitated  droplet  enhances  the  mass  transfer  in  the  liquid  phase  to  a great  extent. 

In  summary,  if  liquid  phase  mass  transport  is  rate  controlling,  the  overall 
reaction  will  be 

(1)  first  order  with  respect  to  [C«]  in  the  solution  reaction  of  diatomic 


(2)  affected  by  the  change  in 


tndition  of  levitated  droplet. 
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However,  the  liquid  phase  transport  is  rarely  found  to  be  a rate  controlling  step  in 
the  gas  - levitated  liquid  metal  reaction  under  electromagnetic  field  because  of  great 
enhancement  of  mass  transfer  due  to  vigorous  stirring. 

Interface  Chemical  Reaction 

In  dealing  with  the  kinetics  of  the  gas/metal  reactions,  interfacial  chemical 
reactions  again  con  be  divided  into  some  substeps.  When  the  forward  and  the 
backward  reactions  in  this  stage  are  accomplished  through  the  same  substeps  but  in 
reverse  sequence,  it  is  called  a reversible  reaction.  On  the  contrary,  in  some  cases, 
the  forward  and  backward  reactions  undergo  the  different  reaction  routes.  This  type 
of  reaction  is  a irreversible  reaction. 


Under  cenain  conditions  of  pressure  and  temperature,  the  reactions  of  metals 
with  diatomic  gases  such  as  H„  N„  and  O,  leads  to  a solution  reaction  which  can 
be  represented  os 

A,(g)  = 2A(diss)  . (2.69) 

For  some  metal-gas  systems,  the  range  of  gas  pressure  in  which  this  equilibrium  is 
established  is  extremely  low  since,  at  high  pressures,  compounds  such  as  nitrides 
and  oxides,  are  formed.  However,  in  systems  like  Fe-N,  Ta-N,  and  Nb-N,  the 
reaclion(2,69)  can  easily  be  stated  over  a wide  range  of  pressures.  Considering  the 
absorption  of  the  gas  A,  in  metal,  the  following  partial  steps  can  be  distinguished  in 
the  inteifacial  chemical  reaction:  (a)  adsorption  of  gas  molecule  on  the  metal  surface 
(physical  adsorption),  (b)  dissociation  of  the  gas  molecules  and  strong  adsorption  of 
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A,(g)  = A, (ad) 

A, (ad)  = 2 A(ad)  (2.70) 

A(ad)  = A(diss)  . 

Dissociative  chemisorption.  Horz  and  Slcinheil[71Horl,  71Hoi2],  Fcng[70Fen] 
and  Rosncr  and  et  al(73Ros,  75Ros]  have  studied  the  nitrogen  absorption  in  solid  Ta 
and  Nbfthin  specimen  with  thickness  less  than  0.5  mm.  where  diffusion  plays  no 
important  role)  at  temperatures  above  1000  °C  in  pressures  ranging  from  10*  to  10* 
Torr.  The  initial  absorption  rate.  v„,  was  found  to  be  proportional  to  the  nitrogen 
pressure  and  dependent  on  reaction  temperature  with  an  activation  energy  of  about 
67  kJ/  mole  N..  The  rate  equation  was  given  as 

r = k P»-  k'CV  = k P„(  1 - ty/C*.1)  (2.71) 

where  k and  k*  are  the  rate  constants  for  the  forward  and  backward  reactions, 
respectively  and  C*.  denotes  the  equilibrium  value  of  the  nitrogen  concentration  in 
the  metal.  These  results  show  that  the  absorption  of  nitrogen  by  Ta  and  Nb 
specimens  is  controlled  by  the  dissociative  chemisorption  of  nitrogen. 

The  kinetics  of  nitrogen  desorption  from  Nb  were  investigated  in  high 
vacuum  for  thin  spccimen(5  < 0.5  mm)  at  temperatures  between  1700  - 2200  °C 

with  high  activation  energies  of  about  519  kl/mole  N. 


indicating  that  recombination  of  nitrogen  atoms  to  molecules  is  the  rate  controlling 


However,  Rosner  and  Coworkers|75Ros]  found  that  dissociation  and 
recombination  are  rate  controlling  for  nitrogen  absorption  and  desorption  of  Ta  only 
in  a definite  range  of  lower  temperatures  and  low  pressures  while  at  high 
temperatures  and  pressures,  either  the  transfer  of  atoms  through  the  surface  or  bulk 
diffusion  becomes  rate  controlling. 


i.  Jehn|70Jeh),  Fromm|7IFrol), 
and  Evans|69Eva]  investigated  the  nitrogen  desorption  in  Mo-N  and  W-N  systems. 
It  was  reported  that  desorption  of  molecular  nitrogen  began  at  a temperature  of 
about  1200  °C  and  the  rate  laws  were  of  the  first  order,  i.e.,  v = - k C„,  with 
relatively  small  activation  energy  of  about  92  kj/mole  N.  Their  analysis  of  the 
kinetic  data  indicated  that  at  the  onset  of  the  reaction,  transfer  of  atoms  from  the 
dissolved  to  the  chemisorbed  state  is  the  rate  controlling  step,  whereas  for  longer 
reaction  times  diffusion  of  the  nitrogen  atoms  to  the  metal  surface  should  be  the 
rate  determining  step  even  with  the  thin  specimens.  It  may  be  concluded  from  the 
nitrogen  desorption  studies  that  the  nitrogen  absorption  in  Mo  and  % 
to  (P*)"1  with  higher  activation  energies  than  found  for  the  Vb  metals.  The  ra 


v = k (P„)'»  - k Q,  = k (P^Hl  - CVQJ  (2.72) 


is  attributed  to  the  rate  control  by  the  transfer  of  nitrogen  atoms  through  the  metal 


surface. 


tbsorptio 


desorption  for  Ta  and 
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other  side  can  be  correlated  with 


dissolution  of  nitrogen  is  highly 
endothermic  for  Mo  and  W. 


Irreversible  reaction 


In  contrast  to  the  metal  - hydrogen  or  nitrogen  systems,  absorption  and 
desorption  of  oxygen  in  many  metals  are  irreversible.  While  oxygen  absorption 
reaction  passes  through  the  same  sequence  of  teacdon  steps  as  discussed  for 


reversible  reactions  earlier,  oxygen  dcsorpuon  is  done  by  the  formation  of  a metal 


xide  into  the  gas 


the  reaction  system.  Considering  only  the  evaporation  of  one  oxide  MeO,  the 
mechanism  of  irreversible  absorption  and  desorption  of  oxygen  can  be  described  by 
the  following  reaction  sequences. 

(a)  0,(g)  = 0,(ad) 

(b)  Oj(ad)  ->  2 O(ad)  (2.73) 

(c)  O(ad)  + Me  ->  MeO(g) 

(d)  0(ad)  = O(diss)  . 


together  with  other  refractory  metal/oxygen  sysla 


ated  to  Nb-O  systen 


At  temperatures  below  1400  - 1600  °C  and  pressures  below  10’  Tore  O,,  the 
reaction  of  Vb  metals,  Nb,  and  Ta,  with  oxygen  results  in  the  solution  of  oxygen 
in  the  metals.  Under  these  conditions,  and  using  thin  specimens,  the  kinetics  of 
absorption  of  oxygen  by  Vb  metals  have  been  investigated  by  many  researchers 
[68Horl,  72Horl,  76Horl,  7SFen,  76Jeh,  65Pas,  65Kof,  88Sch,  77Gra],  For 
sufficiently  short  reaction  times,  the  oxygen  absorption  was  reported  to  be  described 
by  the  following  rate  law|68Horl,  72Froll. 

1+ki,  exp(Q/RT)  . <2  74) 

scheme  described  in  equation  (2,73).  The  activation  energy  of  oxygen  absorption  is 
Q = 47.3  kJ/molcO,  for  V,  59.6  kJ/molcO,  for  Nb  and  73.7  kJ/moloOj  for  Ta. 
Values  of  oxygen  absorption  rate,  v,  decreases  in  the  sequence  V,  Nb,  and  Ta.  The 
rate  cquation(2.74)  implies  that  both  the  transport  of  the  gas  molecules  to  the  metal 
surface  and  the  dissociative  chemisorption  of  the  gas  molecules  arc  rate  determining. 

The  kinetics  of  oxygen  desorption  by  the  evaporation  of  oxides  in  high 
vacuum  have  been  investigated  in  thin  specimens  of  Nb  at  temperatures  above  1400 
°C164Geb,  65Pas,  73Hor].  Desorption  begins  at  temperature  of  1600  °C  for  Nb  with 
the  formation  of  volatile  oxides,  NbO  and  NbO,.  At  oxygen  concentrations  of  less 
than  1 at.%0,  the  evaporation  of  NbO  predominates  while  at  high  oxygen 
concentrations,  the  evaporation  rates  of  NbO  and  NbO,  are  comparable,  leading  to 
the  rale  law  as 

v = k,  C.  + k,  C.1 


(275) 


rate.  G,  Horz 


a!I70Hor,  68Horl,  68Hor2]  found  in  Nb-oxygen 


reaction  coefficient,  r,  which  is  the  fraction  of  the  amount  of  oxygen  dissolved  in 
Nb  among  the  oxygens  delivered  from  the  bulk  gas  to  the  interface,  is  related  to  the 
surface  coverage  by  the  following  expression. 

r = 1 + 1C  exp  (Q/RT)/(l-0)‘  (2-76) 

where  a is  a constant  near  unity.  From  equation  (2.76),  the  data  on  the  surface 
segregation  of  oxygen  on  the  Nb  surface  can  be  obtained  because  the  initial  value 
of  r.  is  available  for  short  reaction  times  with  8 « l[77Gra.  77Kla).  Such  an 
evaluation  of  the  kinetic  data  is  based  on  the  facts  that,  during  oxygen  absorption, 
quasi-equilibrium  between  the  dissolved  and  the  adsorbed  is  established.  For  the  Nb- 
O system,  the  value  of  segregation  enthalpy  was  obtained  to  be  - 50  to  - 63 
fcj/molc  O [77Gra,  77Kla], 

The  reaction  coefficients  for  nitrogen  - metal  and  oxygen  - metal  systems 
arc  plotted  loganthmaiically  as  a function  of  1/T  in  Fig.  2.18.  It  can  be  seen  that  in 
the  oxygen  systems,  the  values  of  reaction  coefficient  becomes  smaller  in  the 
sequence  V,  Nb,  and  Ta.  The  reaction  coefficient  is  much  larger  in  the  oxygen- 
niobium  system  than  in  the  nitrogen-niobium  system,  indicating  that  niobium-oxygen 
solution  reaction  takes  place  much  faster  than  niobium-nitrogen  reaction. 


K.  Ono  and  et  al.  [80Ono]  investigated  the  deoxidation  and  dccarburization 
of  Nb  through  CO  evolution  for  the  samples  containing  various  combination  of  C 
and  O contents  at  temperatures  between  2073  and  2473  K and  at  the  CO  partial 
pressure  of  1 .3*  10‘  Pa.  Their  results  showing  the  carbon  and  oxygen  contents  before 


Kg.  2.18  Reaction  coefficients  for  the  absorption  of  oxygen  or  nitrogen  in 
the  system  V-O,  Nb-O,  or  Ta-O,  or  Nb-N  as  a function  of  the 
reciprocal  absolute  temperature!  thin  samples,  solid  solution 
range,  low  gas  partial  pressures,  sufficiently  short  reaction 
times).  (72HorlJ 


between  2- 10*  and  5-10’  torrl72Horl,  72Hor4,  73Hor2], 
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Fig.  2.19  Correlation  of  the  carbon  and  oxygen  levels  in  niobium  before 
and  after  the  vacuum  degassing  in  solid  statc.[80Ono| 
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Fig.  2.20  Correlation  of  the  carbon  and  oxygen  levels  in  niobium  before 
and  after  melting  in  the  electron  beam  furnace  for  600  seconds. 
[80Ono] 
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Oxygen  absorption 


Decarburization 


IV  Delivery  0/  carbon 


Fig.  2.21  Reaction  scheme  for  ihe  decarburization  of  Vb  mentis  in  flowing 
oxygen. |72Hor] 


CHAPTER  3. 
EXPERIMENTAL 


All  experiments  in  this  study  were  performed  by  electromagnetic  levitation 
technique.  This  technique  is  particularly  suitable  for  studying  gas-Nb  melt 
interactions  because  niobium  with  high  melting  point  and  its  reactive  characteristics 
can  be  melted  without  contact  with  a crucible,  preventing  contamination  by  crucible 
materials.  In  the  section  to  follow,  the  principles  of  this  technique  will  be  discussed 


When  a conductor  is  placed  in  an  alternative  magnetic  field,  it  experiences 
heating  and  force  due  to  the  induction  of  eddy  currents  inside  it.  Sample  heating 
results  from  eddy  current  loss  while  force  arises  from  Lorcntz  magnetic  forces  due 
to  the  interaction  between  eddy  currents  in  the  sample  and  applied  currents  in  the 
coil.  This  force  acts  in  such  a direction  so  as  to  have  the  conductor  move  from  a 
stronger  to  a weaker  part  of  the  magnetic  field.  So,  when  a metal  is  heated  in  two 
coaxial  coil  with  conical  shape  as  shown  in  Fig.  3.1,  a levitation  force  is  generated 
by  the  field  strength  of  bottom  coil  decreasing  in  the  upward  vertical  direction 
[65Fro3,  65Pei,  71Boc,  84Eth|.  At  the  same  time,  a lateral  stability  results  from  the 
field  strength  which  decreases  radially  toward  the  field  axis  and,  in  turn,  produces  a 
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Fig.  3.1  Schematic  diagram  of  levitation  coil 
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restoring  force  toward  the  center  of  the  coil.  The  upper  oppositely  wound  coil 
together  with  a gravitational  force  keeps  the  sample  forced  down  into  the  bottom 
coil,  making  the  sample  heated  with  efficiency  and  stability. 

Electromagnetic  levitation  technique  has  the  following  advantages  for 
studying  gas-metal  reaction  at  high  temperature  [60Bcg,  65Froll. 

1)  The  specimen  is  melted  while  suspended  in  the  space  thus  eliminating 
contamination  from  a crucible  or  container. 

2)  Good  homogeneity  is  achieved  in  the  melt  because  of  the  vigorous  stirring 

3)  No  considerable  temperature  difference  exist  in  the  specimen. 

4)  Melting  is  generally  very  rapid,  often  occurring  in  less  than  a minute. 


Materials 

Niobium  used  in  this  study  was  obtained  from  Cabot  Corporation  in  the  form 
of  rod  with  0.635  Cm  diameter.  The  composition  analysis  is  given  in  Table  7. 
Niobium  rods  were  sheared  into  cylindrical  pieces  weighing  2.1  +0.1  gram.  Close 
control  of  specimen  weight  ensured  that  the  nearly  same  surface  area  existed  in  all 
the  experiments.  The  surface  of  each  specimen  was  cleaned  with  acetone  to  remove 
contaminants  like  grease  and  oil  which  might  be  introduced  inadvertently  during 
weighing  and  handling. 

The  purities  of  gases  used  in  this  investigation  are  shown  in  Table  8.  All 
gases  were  dried  by  passing  them  through  drierite  columns.  Trace  amounts  of 
oxygen  in  the  inert  gases  were 


eliminated  by  passage  through  a column  filled  with 


Tabic  7.  Chemical  analysis  of  Nb, 


Element  N O C 


Fe  Nb 


Concentration  15-30  85-125  15-25  5-10  395-790  10-15  0-10  10-25  Bal. 


Tabic  8.  Composition  analysis  of  gases 


O.  CO,  H,0 


THC  indicates  total  hydrocarbon. 
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titanium  sponge  held  at  800  DC.  The  CO,  impurity  present  in  CO  gas  was  removed  by 
passing  the  gases  through  ascarile  columns. 


Apparatus 

Fig.  3.2  shows  the  details  of  the  experimental  set-up  used  in  this  study.  The 
power  source  for  levitation  was  Lepel  10  kW  R.F.  generator.  It  was  normally 
operated  at  25  - 60  percent  of  its  maximum  power.  The  levitation  coil  was  formed 
&om  1/8  " O.D.,  thin  walled  copper  tubing  wrapped  with  glass  fiber  tape  for 
insulation.  It  had  a conical  shape  with  16  mm  I.D.  at  the  bottom  and  the  full 
conical  angle  of  30  degree.  The  coil  had  14  turns,  the  top  one  opposite  in  direction 
to  the  lower  13  turns  among  which  the  outside  6 turns  are  wound  on  the  inside  7 
turns  in  the  same  direction.  The  reaction  chamber  consists  of  a 5/8  ” O.  D.glass 
tube  open  at  both  ends  and  connected  to  a gas  delivery  system. 

Temperature  was  measured  by  a single  color  pyrometer  calibrated  against  the 
melting  temperature  of  pure  Nb.  The  error  in  temperature  reading  was  + 10  °C.  Some 
degree  of  temperature  control  can  be  obtained  by  adjusting  the  output  power  level. 
In  order  to  increase  the  heating  effects  thereby  increasing  the  sample  temperature, 
the  power  must  be  lowered  to  have  the  sample  rest  lower  in  the  bottom  coil  in 
order  for  more  dense  magnetic  flux  to  go  through  the  sample. 

The  gas  flow  was  measured  using  calibrated  spherical  float  flowmeters  within 
the  accuracy  of  2*.  The  gas  mixture  with  a desired  partial  pressure  of  N,  or  CO 
was  prepared  by  adjusting  the  flow  rates  of  N,  or  CO  and  Ar  gases  at  room 
temperature  using  the  regulating  valves  of  the  flowmeters.  The  rapid  switch  of  gas 
atmosphere  from  the  inert  to  the  reactive  or  vice  versa,  which  is  essential  for 


Fig.  3.2 
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Ihe  prevention  of  time  delay  of  the  reaction  start,  was  done  by  placing  a four-way 
stop  cock  near  the  gas  inlet  to  levitation  column. 


Procedures 

At  the  beginning  of  each  run,  a specimen  is  lowered  to  the  levitation 
position  by  a specimen  positioner.  A flow  of  Ar  and  He  gases  then  was  maintained 
through  the  reaction  column  for  approximately  30  seconds.  After  the  reaction 
column  was  purged  by  inert  gas  mixture,  the  power  was  turned  on,  levitating  the 
specimen.  At  that  time,  the  specimen  positioner  was  raised  up  to  prevent  contact 
with  the  sample  during  melting.  Then  a He  flow  keeps  decreasing  slowly  in  order  to 
increase  the  temperature  of  specimen  up  to  melting  point  of  niobium.  After  about 
5 - 7 seconds  arrest  at  melting  temperature,  the  sample  temperature  keeps  increasing 
and,  finally,  is  stabilized  at  a certain  level  which  depends  upon  the  levitation 
position  of  sample  in  the  coil.  At  this  point,  reaction  temperature  could  be  adjusted 
to  the  desired  temperature  by  controlling  the  power,  thereby  controlling  the 
levitation  position.  Once  sample  temperature  reached  a desired  level,  the  argon  gas 
was  replaced  with  a N,  + Ar  gas  mixture  with  a predetermined  mixing  ratio  almost 
instaneously.  This  was  taken  as  time  zero  for  measuring  the  reaction  time.  After 
reacting  the  sample  with  nitrogen  in  levitated  state  for  a desired  length  of  lime,  the 
sample  was  dropped  and  quenched  on  a copper  chill  plate  by  cutting  off  the  power 
to  the  levitation  coil.  Simultaneously,  the  flow  of  Ar  + N,  gas  was  replaced  with  a 
Ar  + He  flow.  The  total  flow  rate  was  maintained  constant  all  the  time.  This 
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ensured  to  minimize  the  change  in  flow  pattern  in  the  gas  phases  when  nitrogen  gas 
was  introduced. 

The  partial  pressure  of  nitrogen  gas  varied  depending  upon  the  experimental 
conditions.  For  the  study  of  nitrogen  absorption,  the  desired  nitrogen  partial  pressure 
was  obtained  by  adjusting  the  flow  rates  of  N,  and  Ar  gases.  However,  in  studying 
the  effect  of  gas  flow  rate  on  the  absorption  behavior,  a Ar  + 5.14  %N,  gas  mixture 
prcmixcd  in  the  factory  was  used  to  prevent  the  possible  small  fluctuation  of  partial 
pressure  in  each  ran.  Using  the  procedure  described  above,  the  absorption  behavior 
of  nitrogen  in  liquid  Nb  was  investigated  under  various  conditions  of  total  gas  flow 
rate,  partial  pressure,  temperature,  and  time. 

In  studying  the  desorption  behavior  of  nitrogen  in  liquid  Nb,  the  specimen 
was  initially  saturated  with  nitrogen  by  exposing  it  to  N,  + Ar  mixture  with 
nitrogen  partial  pressure  of  0.15  atm  for  4 minutes,  which  was  found  to  be  long 
enough  to  saturate  nitrogen  in  most  of  absorption  conditions  from  the  absorption 
experiments.  The  gas  atmosphere  was  then  changed  to  a pure  Ar  or  Ar+Nj  gas 
mixture  with  lower  partial  pressures  of  nitrogen  instaneously.  Through  such 
procedures,  the  dcsorpdon  kinetics  were  studied  as  a function  of  partial  pressure, 
and  processing  time. 

The  carbon  and  oxygen  absorption/desorption  experiments  were  carried  out  in 
a flow  of  CO  + Ar  gas  mixture  with  various  partial  pressures  of  CO.  In  the  study 
on  the  effect  of  gas  flow  rate,  Ar  + 3%  CO  premixed  gas  was  used.  Most  of 
experimental  procedures  was  basically  the  same  as  those  in  nitrogen  experiments. 
The  carbon  and  oxygen  absorption  kinetics  were  examined  as  a function  of  gas  flow 


CO  pania 


initially  carbon  and  oxygen  were  absorbed  in  Ihe  sample  in  Ihc  CO  + Ar 
atmosphere  with  CO  partial  pressure  of  0.03  atm.  at  2813  K for  4 minutes.  This 

CO  with  lower  partial  pressure  of  CO. 

methods.  Nitrogen  analysis  was  done  (l)by  measuring  the  weight  gain/loss  using  the 
balance  with  the  precision  of  0.001  gram,  and  (2)  by  a flame  atomic  absorption 
method.  Two  methods  gave  essentially  the  same  results  except  at  the  concentration 
higher  than  about  3.5  wt%.  Since  the  composition  data  measured  by  the  flame 
atomic  absorption  method  showed  a scattering  behavior  for  the  duplicate  samples  in 
high  concentration  range,  the  % weight  gain/loss  data  have  been  used  to  analyze  the 

pressure  of  Nb  at  its  melting  point  is  - Iff’  Totr[80Oha|,  the  weight  loss  due  to 

gram,  which  is  negligible  amount  compared  with  the  total  change  in  weight. 
Therefore,  the  entire  weight  gain/loss  could  be  attributed  to  nitrogen  gained/lost. 

On  the  other  hand,  carbon  and  oxygen  contents  of  the  samples  were  analyzed 
by  the  flame  atomic  absorption  technique.  The  detectability  limit  of  this  technique 
was  0.(1001  %. 
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Quadrupolc.  Mass  .Stcnromcig 

Change  of  partial  pressures  of  reactant  and  product  gases  during  the 
processing  was  monitored  by  a Quadrupole  residual  gas  analyzer  manufactured  by 
AMETEKfmodel  MA  200).  Since  the  pressure  of  reacdon  column  was  estimated  to 
be  about  1 atmosphere,  a pressure  reduction  system  was  attached  to  reduce  the 
pressure  to  below  10**  lore,  allowable  pressure  range  by  a Quadupolc  analyzer,  with 
a quartz  capillary  tube  of  1 pm  inside  diameter.  As  soon  as  CO  gas  is  introduced 
into  the  reacdon  column,  monitoring  of  partial  pressure  is  started  by  opening  the 
valve  to  the  analyzing  chamber,  which  was  already  evacuated  to  -10*  Tore  by 
turbomolecular  pump.  The  precision  of  partial  pressure  reading  is  +2  % according  to 
the  manufacturer's  specification.  Using  these  procedures,  cracking  patterns  before 
and  during  the  processing  were  obtained  initially  to  identify  which  gases  are 
involved  in  the  reaction  as  reactant  or  product  gases.  In  the  next  run,  the  change  of 
partial  pressures  of  gases  of  interest  were  monitored  at  a constant  time  interval  to 
examine  the  kinetic  characteristics  of  the  reaction. 


CHAPTER  4 

EXPERIMENTAL  RESULTS 
Nitrogen 

For  the  determination  of  solubility  and  the  mechanism  of  nitrogen  - liquid 
Nb  interaction,  experiments  were  performed  to  study  the  effects  of  partial  pressure 
of  nitrogen,  temperature,  and  time  on  the  nitrogen  absorption  and  desorption 
behavior.  Prior  to  these  experiments,  it  was  examined  how  gas  flow  rate  affects 
absorption  kinetics,  providing  the  idea  of  the  range  of  flow  rate  in  which  the  effects 
of  the  above  parameters  can  be  studied  independently  of  the  gas  flow  rate.  In  the 
following  section,  the  results  will  be  dealt  with  in  detail. 

Effect  of  Gas  Flow  Rate  on  Nitrogen  Absorption 

The  mass  transfer  rate  in  the  gas  boundary  layer  surrounding  the  melt  is 
influenced  by  the  flow  rale  of  the  gas  stream  as  described  in  the  previous  chapter. 
In  order  to  examine  the  effect  of  the  gas  flow  rate  on  the  nitrogen  absorption 
kinetics,  experiments  were  carried  out  at  various  flow  rates  of  Ar  + N,  mixture 
ranging  from  500  ml/min  to  3800  ml/min,  (60-456  in  Reynolds  number),  all  of 
which  belong  to  laminar  flow  regime  in  the  present  experimental  set-up.  The 
nitrogen  absorption  data  at  the  fixed  processing  times  arc  plotted  as  a function  of 
the  gas  flow  rate  in  Fig.  4.1.  In  all  cases,  the  nitrogen  partial  pressure  and 
temperature  were  fixed  at  0.0514  atm  (5.21-10'  Pa)  and  at  2843  K.  It  can  be  seen 
that  the  amount  of  nitrogen  absorbed  increases  with  increasing  the  flow  rate  in 


Fig.  4.1  Nilrogen  absorption  behavior  of  liquid  Nb  as  a function  of  gas  flow 
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lower  flow  regime  ai  each  fixed  processing  time.  However,  at  the  flow  rates  higher 
than  2000  ml/min,  the  change  of  nitrogen  concentration  becomes  negligible  with 
further  increase  of  the  gas  flow  rate.  Such  behavior  supports  the  view  that  under  the 
present  experimental  conditions  the  gas  phase  resistance  to  the  overall  nitrogen 
absorption  process  becomes  insignificant  at  the  flow  rates  higher  than  2000  ml/min. 


The  effect  of  reaction  temperature  on  absorption  kinetics  was  investigated  for 
three  temperatures,  2743,  2843,  and  2943  K.  At  each  temperature,  the  absorption 
data  were  obtained  for  various  lengths  of  time  under  a constant  flow  rate  of  4000 
ml/min  of  Ar  + N,  mixture  but  with  different  nitrogen  partial  pressures  ranging 
from  0.0125  to  0.25  atin.<1.27-10’  Pa  to  2.54-I0*  Pa).  Fig.  4.2  shows  the  variation 
of  nitrogen  concentration  with  time  for  different  processing  temperatures  at  fixed 
nitrogen  partial  pressures  of  0.25,  0.15,  and  0.05  atm.  Nitrogen  absorption  is  rapid 
at  the  initial  stage  but  the  rate  decreases  with  increasing  processing  time  as  the 
sample  becomes  saturated  gradually.  The  saturation  time  was  around  4 minutes  for 
most  samples.  The  saturation  concentration  decreases  with  increasing  processing 
temperature,  indicating  that  the  process  of  nitrogen  absorption  in  liquid  Nb  is 
exothermic,  as  will  be  discussed  later.  This  finding  is  in  good  agreements  with  the 
observation  that  the  sample  temperature  goes  up  sharply  from  the  stabilized 
temperature  above  melting  point  at  the  instance  the  N\  gas  is  introduced  into  the 
levitation  chamber  as  shown  in  Fig.  4.3. 

In  order  to  investigate  the  effect  of  nitrogen  partial  pressure,  the  data 


obtained  at  2843  K under  various  partial  pressures  of  N„  namely,  0.0125,  0.05, 
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Absorption  Time,  sec 

(a) 

Fig.  4.2  Effect  of  temperature  on  nitrogen  absorption  in  liquid  Nb 
(a)  Pa  = 0.25  atm  (b)  P„  = 0.15  atm  (c)  P*  = 0.05  atm 
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Absorption  Time,  sec. 


Absorption  Time, 
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reaclant  gas 


Fig-  4.3  Typical  temperature  profile  during  nitrogen  absorption 
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0.15,  and  0.25  atm,  are  plotted  against  the  processing  time  in  Fig.  4,4.  It  can  be 
seen  that  the  increase  in  nitrogen  partial  pressure  results  in  more  rapid  nitrogen 
absorption  at  the  initial  stage  and  a shoitcr  time  for  nitrogen  saturation.  The 
saturation  level,  being  dependent  on  the  nitrogen  partial  pressure  in  the  gas, 
increases  as  the  partial  pressure  increases.  Similar  trends  were  also  observed  at  two 
other  processing  temperatures,  2743  K and  2943  K. 


Fig.  4.5  shows  the  typical  nitrogen  desorption  behavior  for  three  nitrogen 
partial  pressures  at  the  reaction  temperature  of  2943  K.  In  this  figure,  the  initial 
concentration  is  taken  as  the  saturated  concentration  which  is  obtained  by  levitating 
the  samples  at  2843  K under  the  nitrogen  partial  pressure  of  0.15  atm  for  4 min.  As 
discussed  earlier,  the  gas  flow  was  changed  into  pure  Ar,  0.0125  atm  N.,  or  0.05 
atm  N,  atmospheres  to  make  a desorption  condition  just  after  charging  nitrogen  in 
the  sample.  It  is  to  be  noted  that  the  nitrogen  desorbs  rapidly  at  the  initial  stage  and 
desorption  rate  levels  off  as  degassing  reaction  goes  on.  The  amount  of  nitrogen 
desorbed  at  a given  processing  time  increases  as  nitrogen  partial  pressure  in  the  gas 
stream  decreases. 


Carbon  Monoxide 

In  order  to  determine  the  mechanism  and  the  reaction  rate  of  CO-molicn  Nb 
interactions,  in  addition  to  pressure-temperature-concentration  relationship  in  the 
saturated  state,  the  absorption/desorption  behaviors  of  carbon  and  oxygen  in  CO 
atmosphere  were  studied  under  various  combination  of  the  processing  variables,  i.e. 
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Fig.  4.4  Effect  of  nitrogen  partial  pressure  on  nitrogen  absorption 
in  liquid  Nb 
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Fig.  4.5  Typical  nitrogen  desorption  behavior  in  liquid  Nb 
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partial  pressure  of  CO,  temperature,  gas  flow  rate,  and  the  processing  time.  In  the 
sections  to  follow,  the  effect  of  each  parameter  will  be  presented  in  detail. 

Effect  of  CO  Partial  Pressure  on  C and  O Absorption 

The  effect  of  CO  partial  pressure  on  absorption  behaviors  of  carbon  and 
oxygen  in  liquid  Nb  were  studied  in  Ar  + CO  streams  with  different  CO  partial 
pressures.  In  this  experiment,  temperature  and  flow  rate  were  fixed  at  2813  K and 
4000  ml/min.  Absorption  data  for  carbon  and  oxygen  are  shown  as  a function  of 
processing  time  in  Fig.  4.6  and  Fig.  4.7,  respectively.  Data  point  corresponding  to 
the  same  processing  conditions  in  both  figures  are  taken  from  the  identical  sample. 
It  can  be  seen  that  carbon  and  oxygen  are  dissolved  into  liquid  Nb  together.  This 
behavior  indicates  that  the  reaction,  CO(g)  = (CJ  + [O],  plays  an  important  role  in 
the  dissolution  of  C and  O into  Nb  from  CO  gas.  It  can  also  be  shown  that  carbon 
is  saturated  at  a certain  concentration  depending  upon  the  CO  partial  pressure.  The 
saturation  takes  places  at  a higher  concentration  of  C and  in  a shorter  time  as  the 

minutes.  Similar  behavior  is  also  demonstrated  for  the  oxygen  absorption  but  the 
saturation  level  is  lower  compared  with  the  C saturation  concentration  at  a given 
CO  partial  pressure. 

EffcCLOf  Temperature  on  C and  O Absorption 

Fig.  4.8  and  Fig.  4.9  show  the  absorption  behavior  of  carbon  and  oxygen  at 
different  reaction  temperatures  under  CO  partial  pressure  of  0.125  atm  and  flow  rate 
of  4000  ml/min.  At  all  reaction  temperatures,  liquid  Nb  is  saturated  with  carbon  and 
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Fig.  4.6  Effect  of  CO  paitial  pressure  on  carbon  absorption  in  liquid  Nb 
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Fig.  4.7  Effect  of  CO  partial  pres 


Fig.  4.8  Effect  of  temper 


absorption  from  CO  in  liquid  Nb 


sygen  absorption 


CO  in  liquid  Nb 
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order  to  examine  ihc  interrelationship  between  carbon  and  oxygen  concentrations 
during  absorption,  concentration  data  of  the  samples  used  in  Fig.s  4.6-4.9  are  plotted 
in  Fig.  4.10.  In  this  plot,  C and  O concentrations  marked  by  each  point  were  taken 
from  the  identical  sample.  It  is  of  interest  that  there  exists  quite  a good  linear 
relationship  between  oxygen  and  carbon  concentrations  dissolved  in  liquid  Nb, 
indicating  that  CO(g)  = [C]  + [O]  plays  a major  role  in  the  dissolution  reaction. 

reaction  takes  place  simultaneously  as  a side  reaction,  resulting  in  the  oxygen  loss 

Effect  of  Oas  Flow  Rale  on  C and  Q Absorption 

In  order  to  check  the  effect  of  the  gas  flow  rate  on  the  absorption  kinetics  of 
CO  in  liquid  Nb,  the  absorption  behavior  was  studied  by  varying  the  gas  flow  rate 
in  the  range  of  3000  to  6000  ml/min  (360  to  720  in  Reynolds  number)  with  CO 
partial  pressure  and  temperature  fixed  at  0.03  atm  and  2813  K,  respectively.  Fig. 
4.11  and  Fig,  4.12  summarize  the  variation  of  C and  O concentrations  with  the  gas 
flow  rate  at  fixed  processing  times.  The  amounts  of  C and  O absorbed  are  affected 
slightly  at  lower  flow  rates  in  all  cases  but  at  the  flow  rates  higher  than  4000 

the  gas  flow  rate.  Such  behavior  illustrates  that  the  transfer  rate  of  CO  in  the  gas 

higher  than  4000  ml/min.  Similar  observation  was  made  for  N,  - Nb  absorption 
kinetics.  Otherwise,  absorption  kinetics  would  have  shown  the  dependence  on  the 


CARBON  CONC.,  at  % 


Fig.  4.10  Interrelationship  between  oxygen  and  carbon  concentrations  during 
CO  absorption 
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Flow  Rate(ml/min) 


Fig.  4.11  Carbon  absorption  behavior  of  Nb  from  CO  gas 


function 
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Fig.  4.12  Oxygc 


of  Nb  from  CO  gas 
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gas  flow  rale  since  the  thickness  of  boundary  layer  decreases  with  increasing  ihe 
flow  rate.  Thus,  it  should  be  concluded  ihai  under  the  present  experimental 
conditions  the  gas  phase  transport  is  not  the  rate  controlling  step. 

Gas  Phases  Analysis  by  Quadrupole  Mass-Specirotnclcr 

During  the  CO  absorption,  the  partial  pressures  of  CO  and  CQi  of  outlet  gas 
were  determined  continuously  using  a quadrupole  mass  spectrometer.  Fig.  4.13 
shows  the  change  of  partial  pressures  of  CO  and  COs  normalized  against  die  input 
CO  partial  pressure  with  processing  time.  It  can  be  seen  that  over  the  entire 
processing  period  COj  pressure  remains  nearly  constant  at  about  5%,  the  value  of 
which  is  the  same  as  the  CO,  partial  pressure  of  the  residual  gas  in  the  analyzing 
chamber  of  Quadrupole  mass  spectrometer.  On  the  other  hand,  the  initial  CO 
pressure  detected  at  the  outlet  is  about  half  the  input  CO  pressure  at  the  onset  of 
CO  introduction  into  reaction  section.  As  the  reaction  proceeds  further,  CO  outlet 
pressure  keeps  increasing,  finally  reaching  the  slate  in  which  CO  pressure  remains 
nearly  constant  Considering  that  the  precision  of  pressure  reading  is  +2  %,  die 
discrepancy  between  the  value  of  CO  input  and  pressure  reading  at  the  final  stage 
seems  to  come  from  the  error  of  pressure  readings  or  not  long  enough  time  for 
saturation.  The  observation  described  above  indicates  that  initially  about  half  of  the 
CO  input  is  removed  from  the  gas  phase  as  it  dissolves  in  the  liquid  metal  but  the 
CO  removal  rate  is  decreased  gradually  as  the  process  proceeds  further,  and  it 
finally  becomes  negligible^  the  height  of  the  hatched  area  is  proportional  to  the  rate 
of  CO  dissolution  into  the  sample  or  the  rate  of  CO  removal  from  the  gas  phase). 
This  behavior  is  in  good  agreement  with  initial  rapid  absorption  of  C or  0 in  the 


TIME  (SEC) 


of  outlet  gas  during  C 


Fig.  4.13  Changes  of  partial  pressures  of  CO  and  C02 
and  O absorption  from  CO  gas. 
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interaction  of  CO  with  liquid  Nb,  finally  followed  by  the  saturation  to  certain 
concentrations  of  C and  O as  shown  in  Fig.  4.6  and  Fig,  4.7.  Thus,  it  can  be 
concluded  from  the  present  results  that  C and  O absorption  is  accomplished  through 
CO  dissolution  and  exclude  the  possibility  of  the  reaction,  2 CO  = CO,  + [C],  in 
the  absorption  of  C and  O from  CO  gas. 


Carbon  and  oxygen  desorption  behavior  was  studied  in  atmospheres  of  pure 
Ar  and  CO  + Ar  mixture  with  the  CO  partial  pressure  of  0.0075  atm,  or  0.015  atm 
at  2813  K for  the  samples  which  were  saturated  with  0.03  atm  CO  gas  at  the  same 
temperature.  Concentration  changes  of  carbon  and  oxygen  during  desotption  are 
shown  in  Fig.  4.14  and  Fig.  4.15,  in  which  the  concentrations  at  aero  time 
correspond  to  the  saturated  ones  at  the  processing  conditions  of  CO  charging  prior 
to  desorption  experiment.  There  occurs  a gradual  decline  in  carbon  or  oxygen 
concentration  but  the  rate  slows  down  as  time  goes  on.  It  is  also  apparent  that  the 
carbon  and  oxygen  desorb  more  rapidly  with  lowering  CO  partial  pressure. 

To  examine  the  effect  of  temperature  on  the  desorption  behavior,  the 
desorption  experiments  were  carried  out  in  a stream  of  pure  Ar  at  2753,  2813,  2873, 
and  2933  K for  the  samples  saturated  under  the  CO  partial  pressure  of  0.03  atm  at 
each  temperature.  Since  the  initial  concentrations  are  different  depending  upon  the 
saturation  temperature,  the  raw  concentration  and  the  concentration  normalized 
against  the  initial  concentration  arc  plotted  together  for  carbon  in  Fig.  4.16  and  for 
oxygen  in  Fig.  4.17.  Typical  desorption  behavior  is  observed  at  each  temperature  in 
that  initial  rapid  desorption  is  followed  by  gradual  decrease  of  desorption  rate.  It  is 
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Fig.  4.14  Effect  of  partial  pressure  of  CO  on  carbon  desorption  in  liquid  Nb 
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TIME(SEC) 


Fig.  4.15  Effect  of  partial  pressure  of  CO  on  oxygen  desorption  in 
liquid  Nb 


TIME(SEC.) 


Fig.  4.16  Effcci  of  temperature  on  carbon  desorption  in  liquid  Nb  under 
(a)  raw  data  (b)normalized  concentration  against  initial 


TIME(SEC-) 


Fig.  4.17  Effect  of  temperature  on  oxygen  desorption  in  liquid  Nb  under 


(a)  raw  data  (b)  normalized  concentration  against  initial 


no 

shown  in  ihe  normalized  concentration  plots  that  the  increase  of  temperature 
enhances  the  desorption  rate  over  the  entire  period. 

In  Fig.  4.18,  the  amount  of  oxygen  desorbed  is  plotted  against  that  of  carbon 
desorbed  in  unit  of  atomic  % for  the  samples  used  in  Rg.  4.16  and  Rg.  4.17  in 
order  to  reveal  how  one  quantity  is  related  to  the  other  during  desorption.  In  this 
graph,  each  point  indicates  one  sample  and  each  symbol  represents  the  data  obtained 
at  the  same  temperature.  The  desorbed  amount  of  each  element  is  calculated  by 
subtracting  the  values  of  concentration  at  each  processing  time  from  the  initial 
concentration.  The  slope,  which  is  indicative  of  atomic  ratio  of  oxygen  to  carbon  in 
desorbed  amounts,  is  approximately  unity  in  all  temperatures.  This  ratio  implies  that 
only  the  reaction  (Cl  + (01  ->  CO(g)  is  responsible  for  the  observed  desorption 
behavior,  excluding  the  possibility  of  oxygen  escape  through  the  evaporation  of 
niobium  oxide  in  a pure  Ar  atmosphere.  If  both  mechanisms  were  working  together, 
the  amount  of  oxygen  desorbed  would  have  been  larger  than  that  of  carbon 
desorbed  in  a given  desorption  time,  and  would  have  shown  a slope  greater  than 
unity.  As  described  previously,  the  samples  at  the  starting  point  of  desorption 
contain  carbon  in  excess  since  the  initial  concentrations  of  carbon  and  oxygen  are 
the  saturated  values  corresponding  to  CO  dissolution  condition.  So,  the  present 
results  agree  well  with  the  results  reported  by  othcrs[80Ono]  in  that  only  the 
desorption  through  CO  evolution  is  responsible  for  C and  O refining  when  the 
excess  carbon  exists  in  the  sample.  In  this  case,  the  excess  amount  of  carbon  will 


ain  after  oxygen  is  depleted. 


Fig.  4.18  Interrelationship  between  the  desorbed  amounts  of  oxygen  and 
carbon  in  atomic  % during  CO  desorption. 


CHAPTER  5 
DISCUSSION 


From  the  thermodynamic  calculations  of  the  formation  of  NbjN  presented 
earlierfsee  Table  4),  it  can  be  seen  that  the  minimum  nitrogen  partial  pressure  for 
the  formation  of  the  compound  is  0.3  atm.  at  2773  K.  Thus,  excluded  is  the 
possibility  of  nitride  formation  on  the  surface  of  liquid  Nb  droplet  under  the 
experimental  condition  used  in  this  study.  Accordingly,  only  nitrogen  solution 
reaction  is  taken  into  consideration  in  dealing  with  the  solubility  and  the  kinetic 
modelling. 


The  solution  reaction  of  nitrogen  in  Nb  may  be  described  by  the  following 


where  a,  is  the  activity  of  N in  liquid  Nb  at  temperature  T.  P„  is  the  partial 
pressure  of  nitrogen  in  the  gas,  and  AG°  is  the  Gibbs  energy  change  for  the  solution 
reaction.  As  discussed  earlier,  if  the  solution  obeys  Henry's  law  and  hypothetical  1 
wt.%  solution  of  nitrogen  is  taken  as  the  standard  state,  the  activity  term  in 


equation: 


1/2  N,(g)  = N fin  liquid  Nb) 
the  equilibrium  constant,  K,  for  the  reaction  is  given  by 
K = — prmr  = exp  (-AG°/RT) 


(5.1) 


(5.2) 
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Equation  (5.2)  can  be  directly  replaced  by  wt.%.  Consequently.  Equation  (5.2) 
becomes  identical  with  Sicverts  law,  and  the  concentration  of  N in  Nb  can  be 
related  to  the  partial  pressure  of  N,  as  follows: 

C„=  e*p(-£?°-  -^T)  P.^kCDPw”  (5.3) 

where  C„  is  the  concentration  of  nitrogen  in  liquid  Nb  in  wt.%,  k(T)  is  the 
temperature  dependent  Sieverts  constant,  and  4H°  and  AS°  are  the  standard  enthalpy 
and  entropy  of  solution,  respectively. 

Based  upon  Equation(5.3),  the  saturated  nitrogen  concentration  at  each 
nitrogen  partial  pressure  and  reaction  temperature  is  plotted  against  the  square  root 
of  nitrogen  partial  pressure  as  shown  in  Fig.  5.1.  A marked  deviation  from  linearity 
is  observed  at  high  nitrogen  concentrations,  indicating  that  Sieverts  law  does  not 
hold  over  the  entire  range  of  temperatures  and  partial  pressures  investigated. 
Nevertheless,  Sieverts  law  can  be  applied  at  low  nitrogen  concentration  (below 
approximately  2 wt.%).  Taking  Sieverts  constant  at  each  temperature,  k,  which  is 
the  slope  of  the  linear  portion  of  the  plots  in  Fig.  5.1,  the  values  of  log  k are 
plotted  as  a function  of  IfV  in  Fig.  5.2.  From  the  slope  and  intercept  of  this  plot, 
the  standard  enthalpy  and  entropy  of  solution  were  calculated  to  be  - 236.4  + 23.3 
kJ/mol  and  - 65.3  + 8.3  J/K  mol,  respectively.  Using  these  quantities,  the  solubility 
of  nitrogen  in  liquid  Nb  can  be  represented  as  a function  of  temperature  and  partial 
pressure  of  nitrogen  as  follows: 

log  Q,  = . 3.41  + 1/2  log  P„  (<  - 2 wt.%)  (5.4) 

where  C,  is  in  wt.%,  T is  in  K,  and  Psj  is  in  atm.  It  must  be  emphasized  that 
Equation  (5.4)  is  valid  only  in  the  range  where  Sicvcns  law  is  applicable.  This  is 
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concentration  as  a function  of  Ps,1"  at  various 


Fig.  5.1  Plot  of  equilibrium 
temperatures 
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Fig.  5.2  Arrhenius  lypc  plot  of  log  k.  Sieven 
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slightly  different  from  the  solubility  data  at  2800  K reported  by  Rcvyakin  and 
al[76Rev).  For  example,  the  solubility  of  N at  2800  K under  0.03  atm  nitrogen 
partial  pressure  is  calculated  to  be  10.3  at.%  according  to  equation  (5.4)  while  the 
calculated  value  at  the  same  conditions  according  to  their  solubility  equation  was 
7.9  at.  %.  The  discrepancy  may  perhaps  be  due,  at  least  in  part,  to  the  errors  in  the 
measurements  of  temperature,  partial  pressure,  and  compositions. 

Reaction  Sequence 

Since  nitrogen  absorption  in  liquid  Nb  is  a heterogenous  process,  it  can  be 
considered  to  occur  in  the  following  five  disdnet  elementary  steps,  as  shown 
schematically  in  Fig.  5.3.  Each  elementary  step  and  its  rate  equation  can  be 
described  as  follows: 

(1) gas  phase  transport 

(i)transport  of  nitrogen  molecules  from  the  gas  phase  to  the  gas/metal 
Nj(b)  ->  N,(i) 

r = k.  (A/VKPim,  - P*,,,)  (5.5) 

(2) chemical  reactions 

N2(i)  + v <=>  N,* 

r = (A/V)(k,  Pw„  C,  - k,'  C*,)  (5.6) 

(iii)dissociadon  of  the  adsorbed  molecules  into  adsorbed  atoms. 

N,*  + v <=>  2 N* 


= (A/V)(  k,  Qa.  C,  - k,’  <V) 


(5.7) 


117 


Fig.  5.3  Reaction  mechanism  of  nitrogen  absorpuon/dcsorpuon  in  liquid 


I Gas  phase  transport 

II  Adsorption  of  nitrogen  molecule  on  the  liquid  surface 
ill  Dissociation  of  adsorbed  nitrogen  molecules  into  atoms 

IV  Incorporation  of  the  adsorbed  atom  to  subsurface  layer 

V Liquid  phase  transport 
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(iv)incorporation  of  the  adsorbed  atoms  to  subsurface  layers,  and 
N*  <=>  NO)  + v 
r = (A/VXk.  Q,.  - k/C^C.)  (5.8) 

(3)liquid  phase  transport 

N(i)  ->  N(b) 

r = k,(A/V)(C*ra  - Q^)  (5.9) 

surface,  v and  C,  stand  for  a vacant  site  available  for  reaction  and  the  concentration 
of  vacant  site,  respectively.  The  * indicates  the  entities  adsorbed  on  the  surface,  k,, 
k,,  k„  and  k,’  are  the  mass  transfer  coefficient  in  the  gas  phase,  the  mass  transfer 
coefficient  in  the  liquid  phase,  the  rate  constant  for  forward  reaction,  and  the  rate 
constant  for  backward  reaction,  respectively,  A is  the  specimen  surface  area,  and  V 
is  the  melt  volume.  It  is  to  be  noted  that  the  rates  of  elementary  steps  (ii),  (iii)  and 
(iv)  are  expressed  as  the  difference  between  a rate  in  the  forward  reaction  and  an 
opposing  rate  because  each  elementary  reaction  is  reversible.  In  contrast,  rate 

based  on  the  fact  that  both  processes  are  basically  diffusion  controlled  in  the 
boundary  layer. 

the  various  elementary  steps  list  above,  the  slowest  one  or  a combination  of  steps 


will  control  the  overall  reaction  rate.  Thus,  characteristics  of  each  step  will  be 


_Sle 


If  the  rate  of  nitrogen  absorption  is  controlled  by  the  gas  phase  transport, 


concentration  or  directly  proportional  to  the  nitrogen  partial  press 

2.  The  rate  of  nitrogen  absorption  would  be  influenced  by  the  flow 
input  gas  stream  because  gas  phase  mass  transfer  coefficient  is  a 
of  Reynold  number  which  includes  the  velocity  term. 

3.  The  apparent  activation  energy  would  be  small,  i.e.,  between  1 at 
kcal/mol.(68Med] 


If  the  adsorption  of  nitrogen  molecules  on  the  surface,  step(ii),  is 
controlling  step,  then  the  overall  reaction  rate  will  be  dominated  by 


nitrogen  molecules  adsorption.  In  other  words,  all  other  steps  are  vi 
equilibrium  or  very  rapid.  The  rate  equation  in  this  case  can  be  written  as 


r = (A/V)(k,  Pw  C.  - k,'  C„‘) 
By  assumption,  it  implies  that 


(1)  gas  phase  transport  is  very  rapid. 


PfMM  = = P»  (5.10) 

K,  = Me,'-  (CVVCVQ)  (5-11) 

(3)  dissolution  step  is  also  at  virtual  equilibrium. 

K,  = W = (CroCyCs.)  (5.12) 

(4)  liquid  phase  transport  is  very  fast. 

Qm  = Q«„  =C,  (5.13) 

where  ’K,'s  are  equilibrium  constants  for  each  substep.  Substituting  equations 
(5.10)  - (5.13)  into  equation  (5.6),  it  yields 

r = (A/V)(k,  P„  C.  - (kj'/K,  KJC,  C,')  . (5.14) 

If  the  equilibrium  constant  for  overall  nitrogen  absorption  is  denoted  by  K,  then 

K = K,  K,  K,!  = C*.’  / PV|  . (5.15) 

Using  equation  (5.15),  equation(5. 14)  can  be  further  reduced. 

r=  (A/VOtk^OaO,.1  - Q,1)  (5.16) 


Equatton(5.16)  is  the  rate  equation  in  case  the  adsorption  of  nitrogen  molecule  on 
the  surface  is  the  rate  controlling  step.  It  follows  from  the  equation  that  the  rate  of 
nitrogen  absorption  must  exhibit  a second  order  dependence  with  respect  to  nitrogen 
concentration  in  the  melt  if  the  overall  reaction  rate  is  indeed  controlled  by  step  (ii). 


When  the  dissociation  of  adsorbed  nitrogen  molecules  into  atoms  is  the  rate 
controlling  step,  the  rate  of  nitrogen  absorption  is  given  by 


= (A/V)( 


k>’  (V)  . 


(5.7) 
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In  a similar  way,  the  inlerfacial  and  adsorbed  quantities  may  be  eliminated  from 
equation(5.7).  The  rate  equation  is  then  written  as  follows: 

r = ( A/VX  k',C.’/KJ«V  - C„')  (5.17) 

From  the  above  equation,  it  is  seen  that  the  rate  of  nitrogen  absorption  is  also 
second  order  with  respect  to  nitrogen  concentration  if  overall  reaction  is  controlled 


The  adsorbed  nitrogen  atoms  must  dissolve  in  the  liquid  phase  in  order  for 
the  nitrogen  absorption  to  proceed.  If  the  dissolution  of  nitrogen  atoms  is  the  rate 
controlling  step,  then  the  rate  equation  is. 

r = (A/V)(k,  C*.  - k,'  Cg,,,  C.)  (5.8) 

The  elimination  of  the  interfacial  and  adsorbed  quantities  from  equation  (5.8)  by 
using  the  concept  of  virtual  equilibrium  results  in  the  following  equation. 


r = (A/V)  k/C.  (Cs.  - C*)  (5.18) 

According  to  equation  (5.18),  the  overall  reaction  is  fitst  order  with  respect  to 


If  the  nitrogen  absorption  is  controlled  by  the  liquid  phase  transport,  the  rate 
of  nitrogen  absorption  will  then  be  given  by 
r = (A/V)  k,  (C*.  - QJ  . 


(5.9) 
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This  equation  shows  a first  older  dependence  with  respect  to  the  nitrogen 
concentration.  Also,  liquid  phase  mass  transfer  coefficient,  k,.  will  be  affected  by  the 
change  of  hydrodynamic  condition  in  the  liquid  droplet. 

So  far,  several  rate  equations  have  been  developed  for  the  elementary  step 
which  is  assumed  to  be  rate  controlling.  Each  rate  equation  permits  certain 
predictions  which  must  be  realized  if  the  elementary  step  underlying  the  rate 
equation  is  indeed  the  rate  controlling  step.  Those  are  summarized  in  Table  9. 

The  Older  of  Reaction  and  Rate  Controlling  Step 

Based  upon  the  characteristics  of  each  elementary  step  as  a rate  controlling 
step,  it  is  possible  to  determine  the  order  of  reaction  and  the  rate  controlling  step  in 
the  nitrogen  absorption  reaction.  As  can  be  seen  in  Table  9,  the  reaction  will  have 

depending  on  the  elementary  step.  The  two  types  of  process  can  be  represented  as 
follows: 

First-order  reaction: 

^ = 4 (5.19a) 


or  integrated  form  of  equation(5.19a)  is 


= kr  4-  ( V > (5.20a) 

or  integration  of  equation(5.20a)  gives 
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Tabic  9.  Summary  of  ihe  characteristics  of  each  elementary  step  as  a rate 
controlling  step 


Elementary  Flow  rate  nitrogen  nitrogen  hydrodynamics 
step  partial  pres.  cone.  in  liquid 


Gas  phase 

transport  O «=  PK,  C*1  X 


Adsorption  of 

N,  molecules  X « P„  “ Cs2  X 


Dissociation 

x - p*  - <v  x 


Dissolution 

X « p„w  «C,  X 


Liquid  phase 

transport  X •*  Ps,'”  =•  Q,  O 


X;  no  effect 


(5.20b) 


2T^T 


[In  ( 


Qt±Qt_ 


c*-c, 

where  On,  C*,,  and  G.,  are  nitrogen  concentration  at  time  «.  at  t = 0.  and  at 
equilibrium,  respectively.  A is  the  specimen  surface  area,  V is  the  melt  volume, 
and  k,  and  k3  arc  the  first-and  second-order  rate  constant,  respectively.  Now,  if  gas 


te  will  be  second  order  and 
rs  transfer  coefficient,  k0, 


e calculated  and  plotted 
:n  that  the  rate  constant 


phase  transport  (step(i»  is  rate  controlling,  the  overall  re 
be  affected  by  the  change  of  the  gas  flow  rate  because  m 
is  influenced  by  the  gas  velocity.  Accordingly,  with  the  s 
4.1,  the  second-order  rate  constants  for  each  flow  rate  wi 
as  a function  of  Reynolds  number  in  Fig.  5.4.  It  is  s 
increases  with  increasing  gas  flow  rate  up  to  approximately  200  in  Reynolds  number 
beyond  which  it  becomes  essentially  independent  of  the  gas  flow  rate.  This  indicates 
clearly  that  at  lower  flow  rates,  the  reaction  rate  is  influenced  by  the  gas  flow 
velocity,  implying  that  the  gas  phase  mass  transport  controls  the  overall  reaction 
rate.  Physically,  in  this  regime,  gas  boundary  layer  is  so  wide  that  the  concentration 
gradient  in  this  layer  is  quite  small,  leading  to  the  slower  reaction  rate.  As  the  flow 
rate  increases  further,  the  mass  transport  rate  of  nitrogen  molecules  keeps  increasing 
because  the  concentration  gradient  is  increased  further,  finally  reaching  a point  at  = 
200  in  Reynolds  number  where  the  mass  transfer  rate  in  the  gas  phase  is 
comparable  to  that  of  the  other  elementary  steps.  Beyond  that  point,  the  gas  phase 
transport  takes  place  faster  than  the  other  elementary  steps.  Consequently,  the 
reaction  control  is  transferred  from  step  (i)  to  one  or  more  of  the  other  elementary 
steps,  resulting  in  the  observed  non-dependence  of  the  rate  constant  on  the  total 


Rate  Constant, (wt%' 1 sec.' 
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Fig.  5.4  Effect  of  gas  flow  rale  on  the  second  order  rate  constant  for 
nitrogen  absorption  in  liquid  Nb. 
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flow  rate.  It  may  be  mentioned  that  the  absoiption/desorplion  data  described 
previously  were  obtained  at  total  flow  rate  of  480  in  Reynolds  number. 
Consequently,  step  (1)  is  not  expected  to  have  any  controlling  effect  on  the  overall 

In  case  step  (v)  is  a rate  controlling  step,  any  change  in  hydrodynamics  of 
the  liquid  phase  will  have  a considerable  effect  on  the  rate  of  nitrogen  absorption 
and  rate  equation  should  be  first  order.  It  is  likely  that  vigorous  convective  motion 
present  in  the  levitated  liquid  Nb  droplet  under  electromagnetic  field  will  greatly 
increase  the  mass  transfer  rate.  This  is  supported  by  Rao  and  Lee’s 
experiment!  82Lee]  on  the  diffusivity  of  carbon  in  an  iron-carbon  droplet  where  the 
carbon  mass  transfer  in  the  liquid  phase  was  enhanced  by  about  ten  times  under 
electromagnetic  conditions  compared  to  a stationary  condition.  Consequently,  step 
(v)  is  not  expected  to  be  a rate  controlling  step  in  the  experimental  system  used  in 
this  investigation.  This  is  further  supported  by  these  experiments  which  show  that 
nitrogen  absorption  is  not  a first  order. 

In  determining  the  rate  controlling  step  among  other  three  elementary  steps, 
it  is  necessary  to  compare  the  data  obtained  experimentally  with  theoretical  rate 
laws.  As  summarized  in  Table  9,  if  step  (ii)  or  (iii)  is  the  rate  limiting  step,  then 
the  reaction  is  a second-order  process  with  respect  to  nitrogen  concentration,  while 
if  step  (iv)  controls  the  overall  reaction  rate,  then  the  reaction  is  a first-order 
process.  It  should  be  noted  that  the  rate  constants  k,  and  k,  in  the  equations  (5.19a) 
and  (S.20a)  are  constant  at  a given  temperature  because  three  elementary  steps  of 
consideration  are  chemical  processes.  In  order  to  determine  if  the  overall  reaction  is 
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first-order  or  not,  the  absorption  data  collected  at  2743,  2843,  and  2943  K under 
various  partial  pressures  of  nitrogen  are  plotted  in  Fig.  5.5  using  Equation  (5.19b). 
In  this  figure,  the  slopes  indicate  the  rate  constant,  k,.  It  can  be  seen  that  in  each 
case,  the  data  arc  split  into  three  or  four  distinct  lines,  each  representing  a different 
partial  pressure  of  N,  in  the  gas  mixture.  If  the  overall  reaction  is  first-order,  the 
data  obtained  at  a given  temperature  should  fall  on  one  line  irrespective  of  the 
nitrogen  partial  pressure  because  the  rate  constant  is  determined  by  the  reaction 
temperature  only.  Accordingly,  this  finding  suggests  that  the  rate  controlling  step  for 
nitrogen  absorption  in  liquid  Nb  is  not  first-order.  This  implies  that  the  elementary 
step  (iv)  does  not  control  the  overall  reaction. 

In  order  to  examine  the  possibility  of  a second-order  reaction,  the  left  hand 
side  term  in  Equation  (5.20b)  was  plotted  as  a function  of  (A/V)t  for  each  reaction 
temperature.  This  is  shown  in  Fig.  5.6  for  temperatures  of  2743,  2843,  and  2943 
K.  It  can  be  seen  that  within  the  limits  of  experimental  error,  the  data  in  each  case 
fall  on  a single  straight  line  regardless  of  the  N,  partial  pressure  used.  Thus,  it  is 
concluded  that  either  elementary  step  (ii)  or  (iii),  which  arc  both  second-order,  must 
be  controlling  the  overall  reaction  rate. 

The  data  analysis  for  desorption  experiments  is  exactly  analogous  because 
nitrogen  desorption  is  expected  to  occur  in  revetse  sequence.  For  the  desorption  data 
in  Fig.  4.4,  the  first  order  and  the  second  order  plots  were  drawn  in  Fig.  5.7  and 
Fig.  5.8,  respectively.  In  these  figures,  the  initial  concentration,  C*,.  and  the 
equilibrium  concentration,  C,.,  were  taken  from  the  absorption  data  obtained  under 


processing  conditions.  The  first-order  plot  shows  a splitting  similar  to  that 
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Fig.  5.5  Kinetics  of  nitrogen  absorption  into  liquid  Nb  under  various 
partial  pressure  of  nitrogen  correlated  with  first  order  rate 
equation 


(a)  2743  K (b)  2843  K (c)  2943  K 
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(A/V)t,  cm  1 sec. 


<b) 
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(c) 
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to) 

Fig.  3.6  Kinetics  of  nitrogen  absorption  into  liquid  Nb  under  various 
partial  pressure  of  nitrogen  correlated  with  second  order  rate 
equation 

(a)  2743  K (b)  2843  K (c)  2943  K 
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(b) 
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(c) 
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(A/V)t,  cm'1  sec. 


Fig.  5.7 


Kinetics  of  nitrogen  desorption  from  liquid  Nb  under  various 
partial  pressures  of  nitrogen  correlated  with  first  order  rate 
equation 
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Fig,  5.8 


Kinetics  of  nitrogen  desorption  from  liquid  Nb  under  various 
partial  pressures  of  nitrogen  correlated  with  second  order  rate 
equation 
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analysis,  it  can  be  concluded  that  the  desorption  reaction  is  again  controlled  by 


Since  the  nitrogen  absorption  process  was  found  to  be  second  order,  the 
apparent  activation  energy  for  nitrogen  absorption  can  be  deduced  from  the 
temperature  dependence  of  the  second  order  rate  constant.  The  Arrhenius-type  plot 

energy  was  calculated  to  be  327.2  ± 20.6  kJ/mole  from  the  slope  and  preexponential 
constant,  k„,  was  2. 29*10'  from  the  intercept.  Therefore,  the  rate  equation  can  be 
expressed  as 


( “T?  Z29°'Uf  "p  ( ' > V I5-21) 

for  2700  K to  3000  K. 

where  C„  and  Q.  are  the  nitrogen  concentration  and  equilibrium  nitrogen 
concentration,  respectively,  in  wt  %,  A is  the  surface  area  in  Cm1,  V is  the  volumn 
in  Cm1,  T is  the  temperature  in  K. 

Temperature  dependence  of  the  desorption  rate  can  be  derived  by  making 
use  of  the  values  of  the  solution  enthalpy,  - 236.4  kj/mole,  and  the  second  order 
rate  constant  for  the  desorption  process  at  2943  K,  a value  of  which  was  obtained 
to  be  1.820* 10*  from  the  slope  in  Fig.  5.9.  As  mentioned  earlier,  the  nitrogen 
solution  reaction  in  liquid  Nb  is  reversible  and  found  to  be  exothermic.  Thus,  the 
apparent  activation  energy  for  desorption  should  be  larger  than  that  for  the 
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Fig.  5.9  Effect  of  temperature  on  second  order  rate  constant,  k,  of 
nitrogen  absorption 


absorption  by  the  amount  equal  to  solution  enthalpy.  Accordingly,  the  rate  equation 

< t?-  >fc-  l.68-10’exp(-  -*^00.  ) £ «V  - Q,’)  (5.22) 

for  2700  K to  3000  K. 

From  the  difference  in  activation  energies  for  both  processes,  it  can  be 
predicted  that  desorption  process  takes  place  more  slowly  than  absorption  reaction 
under  the  same  amount  of  concentration  deviation  from  the  equilibrium  stale.  Such 
prediction  is  demonstrated  well  in  that  the  calculated  rate  constant  for  absorption  at 
2943  K is  3.62-10’  which  is  larger  than  that  for  desorption  process  at  Ihe  same 
temperature,  1.820*10°. 


Carbon  Monoxide 

Solubilities  of  C and  O under  CO  Atmosphere 

In  older  to  get  the  pressure-temperaiure-equilibrium  concentrations 
relationship  concerning  CO-  liquid  Nb  interaction,  the  values  of  0/0,,  the 
multiplication  of  C and  O equilibrium  concentrations,  obtained  at  2813  K(from  Fig. 
4.6  and  Fig.  4.7)  were  plotted  as  a function  of  P„  in  Fig.  5.10.  In  view  of  the 
thermodynamics  of  dissolution  of  C and  0,  a linear  relationship  should  arise  when 
activities  of  C and  O are  used  instead  of  the  concentrations  of  each  elements.  A 

coefficients  of  both  elements  are  relatively  constant  in  the  range  of  temperature  and 
partial  pressures  of  CO  used. 
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Pco(atm) 


Fig.  5.10  Relationship  between  equilibrium  carbon  and  oxygen 
concentrations  and  partial  pressure  of  CO 
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in  which  the  value  of  Ln  Cc.'C0yP„  obtained  at  various  temperature  (from  Fig.  4.8 
and  Fig.  4.9)  were  plotted  versus  1/T.  From  the  slope  and  the  intercept  in  this 
graph,  the  P-T-C  relationship  in  the  saturated  state  was  deduced  to  be 


Some  observations  which  might  be  useful  in  developing  a kinetic  model  for 
C and  O dissolutions  in  liquid  Nb  from  CO  gas  are  summarized  as  follows: 

1.  Both  carbon  and  oxygen  atoms  from  CO  gas  arc  dissolved  into  liquid  Nb. 

( Fig.  4.6  and  Fig.  4.7) 

2.  The  dissolution  process  proceeds  with  a relatively  constant  ratio  of  O to  C 

of  less  than  unity(  - 0.76)  over  the  entire  ron.(  Fig.  4.10) 

3.  There  exists  the  saturated  concentrations  of  C and  O at  a given  partial 

pressure  of  CO  and  temperature.(Fig.  4.6  and  Fig  4.7) 

unity.(Fig.  4.18) 

Based  upon  the  above  observations,  it  can  be  realized  clearly  that  CO(gas)  = 
C(in  Nb)  + 0(in  Nb)  is  the  major  reaction  in  the  dissolution  of  C and  O into  liquid 
Nb  from  CO  gas.  However,  the  atomic  O/C  ratio  less  than  the  unity  makes  it 
necessary  to  determine  how  C surplus  (or  0 deficiency)  arises  during  the  dissolution 


ln  Cc  + ln  C„.  - ln  P„  = 2.32- 10VT  - 2.35 


(5.23) 


where  P_  is  in 


Cc,  and  Co,  are  in 
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Fig.  5.11  Relationship  between  Ln(Cc.QJPco)  and  1/T 
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reaction.  Ii  can  be  speculated  that  there  exists  some  side  reaction(s)  which  causes 
the  observed  results.  There  are  two  possible  reactions  to  explain  such  results. 

i)  2 CO  <=»  [q  + CO.  (5.24) 

ii)  Nb  + O ->  NbO 

or  Nb  + 20  ->  NbO,  (5.25) 

If  reaction  (i)  is  taken  as  a side  reaction,  it  would  provide  another  source  for 
carbon  dissolution  which  can  cause  carbon  surplus.  Thermodynamic  equilibrium  of 
reaction  (i)  indicates  that  reaction  (i)  at  2773  K can  lakes  place  in  the  forward 
direction  only  when  P_/P„  is  greater  than  - 2.5*  10’  under  the  present  processing 
conditions.  Since  the  CO  gas  used  in  this  experiment(initially  P_/P„  is  6.7-10")  is 
expected  to  have  much  low  P J PM,  ratio  even  after  passage  through  an  ascaritc 
column,  reaction(i)  can  be  discarded  among  the  possibilities  to  explain  the  observed 
behavior.  In  addition,  this  view  is  also  supported  by  the  observation  that  CO, 
pressure  remains  almost  constant  during  the  entire  processing  period  as  shown  in 
Fig.  4.12.  Next,  suppose  that  reaction(ii)  occurs  as  a side  reaction.  Then,  oxygen 
loss  will  result  from  the  evaporation  of  niobium  oxide(s).  leading  to  the  oxygen 
deficiency  even  if  the  same  amount  of  adsorbed  oxygen  and  carbon  are  generated 

oxides  is  known  to  become  significant  for  oxygen  desorption  through  niobium 


above  approximately  1900  K 
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greater  than  Iff1  atm  [72Frol,  80Fro,  SlSchlJ.  reaction(ii)  seems  to  be  the  most 
possible  side  reaction  to  explain  the  0 1C  ratio  less  than  unity. 

For  the  desorption  process,  on  the  other  hand,  the  atomic  ratio  of  O/C  of 
about  unity  in  the  desorbed  amounts  indicates  that  simultaneous  escape  of  C and  O 
occurs  through  CO  evolution,  excluding  the  possibility  of  the  decrease  of  oxygen 
concentration  through  the  path  of  oxide  evaporation.  The  difference  in  mechanisms 
between  absorption  and  desorption  can  be  explained  as  follows:  during  the 
absorption,  CO  molecules  are  adsorbed  on  the  Nb  surface  and  dissociated  into  C 

for  the  adsorbed  oxygen,  there  are  two  competing  paths,  say,  dissolution  into  liquid 
Nb  or  desorption  through  niobium  oxide  evaporation.  Of  both  paths,  it  is  surmised, 
based  on  the  observed  O/C  ratio  of  0.76  in  the  samples,  that  the  dissolution  path  is 
dominant  over  the  oxide  evaporation  by  approximately  4 times.  For  desorption 
process,  on  the  other  hand,  the  oxide  evaporation  and  the  recombination  of  C and  O 
may  be  competing  after  the  C and  O atoms  are  transferred  to  the  surface  from  the 
interior  of  liquid  Nb.  The  atomic  O/C  ratio  of  about  unity  in  desorbed  amount 
suggests  that  recombination  process  is  dominant  over  the  oxide  evaporation  to  a 
great  extent  especially  when  carbon  exists  in  excess  initially.  This  behavior  seems  to 
be  closely  related  to  the  equilibrium  vapor  pressures  of  niobium  oxide  and  CO  in 
Nb-C-O  solution  as  shown  in  Fig.  2.6  in  which  the  vapor  pressure  of  CO  is  two  to 
three  orders  higher  than  that  of  niobium  oxides  over  the  niobium  containing  the 
same  amounts  of  C and  O.  It  is  also  to  be  pointed  out  that  the  observed  desorption 
mechanism  is  consistent  with  Ono  and  et  al’s  works|80Ono)  on  the  degassing  of 
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niobium  containing  carbon  and  oxygen  in  which  only  CO  recombination  occurs 
when  carbon  exists  in  excess  whereas  both  CO  recombination  and  niobium  oxide 
evaporation  take  place  together  when  oxygen  exists  in  excess.  Thus,  it  can  be  said 
that  in  absorption,  CO  dissolution  is  responsible  for  C and  O absorption  together 
with  niobium  oxide  evaporation  whereas  in  desorption  only  CO  recombination  path 
is  taken  for  C and  0 desorption.  Taking  the  discussion  done  so  far  into 
consideration,  the  absorption  and  desorption  paths  can  be  visualized  as  shown  in 
Fig.  5.12. 

In  the  followings,  a kinetic  model  on  the  interaction  of  liquid  Nb  with  CO  is 
made  based  upon  the  reaction  CCKg)  <=>  [C]  + [O]  together  with  a side  reaction  of 
niobium  oxide  evaporation  as  shown  schematically  in  Fig.  5.13.  Each  elementary 
step  involved  in  this  scheme  can  be  described  in  the  following  way. 

I.  Gas  phase  transport 


CO*(g) 


CO'(g) 


(5.26) 


II.  Dissociative  adsorption 


CO(g)  + 


2v 


C*  + O* 


(5.27) 


IU.  Dissolution 


a)  Carbon  dissolution 


C* 


C(i)  + 


(5.28) 
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a:  Equilibrium  line 


Fig.  5.12  Reaction  paths  during  CO  absorptiu 


146 


Fig.  5.13  Mechanism  of  carbon  and  oxygen  dissolution  in  liquid  Nb  from 
CO  gas 

I.  gas  phase  transport 

II.  dissociative  adsorption  of  CO 

III.  a.  carbon  dissolution 

b.  oxygen  dissolution 

c.  niobium  oxide  evaporation 

IV.  liquid  phase  transport 


b)  Oxygen  dissolution 
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O* 


20*  + Nb 

IV.  Liquid  phase  transport 


C(i)  ->  [Cl 

0<i)  ->  [O] 


(5.29) 

(5.30) 

(5.31) 

(5.32) 

(5.33) 


Hence,  C*  and  O*  denote  chemisorbed  C and  O,  respectively,  C(i)  and  O(i)  arc 
dissolved  species  near  subsurface  of  liquid  Nb,  v is  a vacant  surface  site,  and  r,  is 
the  reaction  rate  for  an  elementary  step  i.  i,  b,  and  [ 1 indicate  the  near  interface, 
the  bulk,  and  the  dissolved  It  should  be  noted  that  step  IIl-(c)  is  not  included  in 
reaction  sequence  of  desorption  because  C and  O desotpdon  proceeds  only  through 
CO  gas  evolution  as  discussed  previously.  In  this  model,  bound  oxygen  and  carbon 
are  assumed  to  be  chemically  adsorbed  in  an  energetically  common  state.  In  some 
cases,  physisorption  of  gas  molecules  as  a precursor  to  atomic  chemisorption  is 
taken  into  consideration  but  the  current  scheme  combines  both  steps  into  one.  This 
is  because  the  molecularly  adsorbed  gas  intermediates  are  generally  expected  to  have 
very  short  life  time  at  high  temperatures,  providing  no  real  physical  meaning  as  an 
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independent  substep[71Winl,  80Fra).  Another  thing  to  be  noted  is  that  only  two 
oxide  forms  via  NbO(g)  and  NbO,(g)  are  taken  into  account  and  metal  oxide 
evaporation  is  assumed  to  be  an  irreversible  process.  This  assumption  enn  be 
justified  by  the  fact  that  the  evaporating  oxides  are  condensed  on  the  glass  wall 
continuously,  making  the  backward  reaction  virtually  negligible.  Then  the  elementary 
reaction  rate  can  be  expressed  in  terms  of  elementary  rate  constants  and  relevant 


concentrations  according  to  die  mass  action  law  as  fo!lows:I74Ros, ' 
r,  = k,(A/V)(P„'  -P J) 

r,  = kj(A/V)P„'(l  -6c  -9o)! 

r>'  = kjTA/Vjflc  60 

r„  = k„(A/V)9c 

■V  = ko'(A/V)(l  -6c  -0o)  Q,,, 

r.j  = k„(A/V)0o 

V = k„’(A/V)(l  -0c  -0o)  C,*, 

r««  = WA/V)0o 

is*.  = k>tu,(A/V)0o' 

r-  = UA/VXCc,,,  - Cy 

r-  = UA/VNCo,,,  - QJ 

where  0c  and  9o  are  the  fractional  surface  coverage  of  adsorbed  C 


(5.34) 

(5.35a) 

(5.35b) 

(5.36a) 

(5.36b) 

(5.37a) 

(5.37b) 

(5.38) 

(5.39) 
(5.40a) 
(5.40b) 

O.  A and  V 


the  step  i,  and  Cc  and  Co  are  concentrations  of  C and  O species  in  the  sample. 


d0c/dt 


= (A/V)[  k,  PJ(1  -9c  -0O)1  - k,'0c  0o  - k.,  0c 

+ krt'(l  -0c  -0o)  CqiJ  (5.41) 

d0o/dt  = r,  - r,’  - r„,  + r.,'  - 

= (A/V)l  k,  PJ(1  -0c  -0o)‘  - k,'0c  0o  - k„,  0o 

+ k0,’(l  -0c  -0o)  Q,,,  - k»o  00  - kmo,  0o!  ) . (5.42) 

Steady  stare  assumption  of  intermediates  species  gives  d0o/dt  - 0 and  d0c/dt  =■  0. 

As  mentioned  earlier,  the  gas  flow  rate  docs  not  affect  the  C and  O 
absorption  rates  under  the  present  experimental  condition.  Based  upon  the  same 
explanation  applied  to  Nb-  N kinetics,  it  can  be  said  that  the  transport  in  gas 
phase(step  I)  is  not  the  rale  controlling  step,  implying  that  P„“  = P„‘  = P„.  This  is 
because  the  rate  of  mass  transfer  in  the  gas  phase  is  rapid  enough  compared  with 
the  rates  of  other  steps,  resulting  in  no  concentration  gradient  in  the  boundary  layer. 
The  possibility  of  liquid  phase  transport  as  a rate  controlling  step  may  also  be 
excluded  using  the  same  reasoning  as  in  liquid  Nb-N  interaction.  Thus, 
mathematically,  it  can  be  assumed  reasonably  that  Cq,,  = Q,  and  = C0. 

Suppose  that  dissolution  step  for  carbon  (III-a)  and  for  oxygen  (m-b)  are 
the  rate  controlling  step.  The  rate  equations  for  C and  O absorption  can  be  then 
represented  as 


dCydt  = r„  - r„' 

= (A/V)[  lc*  0c  - ka'(l  -9c  -00)  CJ 


(5.43) 
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d(Vdl  = r„  - r„’ 

= (AA')I  k.,  6o  - K,,’(l  -0c  -00)  C„]  . (5.44) 

The  above  assumption  implies  that  dissociative  adsorption  ratefstep  II)  is  so  fast, 
compared  with  the  net  incorporation  rate  of  C(step  Ill-a)  or  0(step  III-b),  that 
stepfll)  is  in  virtual  equilibrium  state.  So,  the  corresponding  adatom  coverage.  0c 
and  0o,  are  virtually  constant  at  the  coverage  of  quasi-equilibrium,  0c-  and  0o', 
throughout  most  of  the  absorption  period.  Therefore,  carbon  and  oxygen  dissolution 

dCc/dt  = (A/V)[  K,  0c1  - ka’(l  -0c'  -00')  QJ  (5.45) 


dCe/dt  = (A/V)[  k.,  0o-  - kj'(l  -0c'  -0o')  CJ  . (5.46) 

On  the  other  hand,  the  net  dissolution  rates  of  C and  O are  zero  in  the  quasi- 
equilibrium  state,  so  that 

dQ/dt  = (A/V)[ko  0C  - l«,’(l  -0c'  -©o')  CJ  = 0 (5.47) 


dCn/dt  = (A/V)[k.,  0o'  - k„,'(l  -0c'  -0o')  CJ  = 0 
where  Q,  and  Cq,  are  the  concentrations  of  carbon  and  oxygen  in  equilibriu 
respectively.  From  equations  (5.47)  and  (5.48),  it  can  be  derived  that 
(1  -0c'  -0o')  = 


(5.48) 


(5.49) 


Inserting  equ 


(5.49)  into  equations  (5.45)  and  (5.46)  giv 
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dQ/dt  = (A/V)[ka  Be-  - kd-  Cc  (ko  Be 

= (AA01  (Q.-Q)] 

= Kc(A/V)(  Cc  - Cc) 


r in  inlegraied  form. 


Cc-Q 


(5.50a) 

(5.50b) 


dQ/dt  = (A/V)[k.,  Bo'  - k„'  Q (k„.  Bo"  / k.,'QJ] 

= (A/V)|  ( Cc  • Co)] 

= Ko(A/V)(  Cc  • C„)  (5.51a) 

or  in  integrated  form, 

■ ln  ( ■%r-cS c > = Ko  ‘ (5’5lb) 

where  Kc  = k.,  Bc’/Q,  and  Ko  = k„,  Bo'/Cc- 

These  equations  (5.50a)  and  (5.51a)  indicate  that  if  step  III  is  the  rate  controlling 
step,  the  overall  reaction  either  for  carbon  dissolution  or  for  oxygen  dissolution  will 
be  first  order  with  respect  to  C or  O concentration,  respectively. 

On  the  other  hand,  let's  assume  that  dissociative  adsorption(step  II)  controls 
the  overall  reaction.  Then  rate  equation  is  written  as 

dQ/dt  = dQ/dt  = r,  - r/  (5.52) 

Hence,  since  dQ/dt  - dQ/dt  and  dQ/dt  + dQ*/dt  + dQWdt  - 
dQ/dt,  the  rate  equations  for  carbon  and 
represented  as 
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dCj/dt  = (A/V)[  k,  P„'(l  -0c  -00)'  - k,'0c  0o]  (5.53) 

dCVdl  = (A/V)[  k,PJ(l  -0c  -0o)’  - kj'0c0o  - (k^  6o  + kh«,0o')]  (5.54) 

The  assumption  of  step  II  being  the  rate  controlling  step  implies  that  dissolurion 
renctionsfstcp  IH-a  and  III-b)  are  so  fast  that  the  virtual  equilibrium  between  the 
adsorbed  and  the  dissolved  species  are  maintained.  Thus,  the  following  relationships 
should  be  satisfied. 

K.,  = k,A„'  = Ccd  - 9c  - 0oV0c  (5.55) 

K„,  = Me.,’  = Co(l  - 0c  - 0o)/0o  . (5.56) 

Substitution  of  the  equations  (5.55)  and  (5.56)  into  the  rate  equation  (5.53)  yields 

dQ/dt  = (A/V)Ikj  PJ1-0C-0O)1  - kj'  Q Cdl-Oc-Ool'/K,.,  KJ  . (5.57) 

Since  the  equilibrium  constant  for  the  overall  reaction  is  related  to  equilibrium 
constants  for  each  elementary  step,  equilibrium  concentrations  of  C and  O,  and 
partial  pressure  of  CO  gas,  we  have 

K = K>  Ka  K.,  = (Mt,')(ka/kc,')(k.1/k<,')  = C*  CJ  P„  . (5.58) 

Substituting  the  relationship  (5.58)  in  equation  (5.57)  gives 

dC/dt  = (A/V)[kj  (l-0c-0o)’(Cc  CoJ/K  - Ic'lWCe  C„)(l-0c-0o)!/K] 

= (kjfKXAAOd-ec-SoWCc,  C„.  - Cc  Co) 

= k,’(A/V)(Cc,  Co.  - Cc  CJ  (5.59) 

where  k,'  is  (k^K)(l-0c-0o)’. 


O/C  atomic 


_ [k,P„‘(l  -6c  -60)*  - lc,'8c9ol  - (k^  0o  + k^Bo1) 

k,  Pj(l  -6c  -00)'  - ki’0c  6o 

(k«»  eo  + k^»8o') 

1 ‘ k,  P„  (!  -6c  -So)’  - k,'0c  So 

= consinnt . (5.60) 

From  die  equation(5.60),  it  can  be  realized  lhal  ( ks„,  0o  + ks„„  So1)  should  be 
always  a constant  fraction  of  [k,  PJ(I  -6c  -6o)!  - ki'Oc  6o].  Thus,  etiuation  (5.54) 


dQ/dt  = K"(AA')(k!  PJ(I  -9c  -0o)'  - k/0c  9o)  (5.61) 

where  K”  is  the  unity  less  the  value  of  the  ratio  (k**,  9o  + kmot6o1)/[k1P„i(l  -0c  - 
0o)!  - k,'0c0o]. 

Thus,  it  is  seen  that  the  rate  equation  for  oxygen  dissolution  is  exactly  the  same 
foim  as  that  for  carbon  dissolution.  Using  the  same  procedure  for  dQ/dt.  equation 
(5.61)  yields 

dQ/dt  = k/(A/V)  ( Cc.  C„.  - Cc  Q, ) (5.62) 

where  k„’  = K”k,’.  Based  upon  the  equations  (5.59)  and  (5.62),  it  can  be  said  that  if 
the  dissociative  adsorption  (step  II)  is  controlling  the  overall  reaction,  both  carbon 
and  oxygen  dissolution  rate  should  be  second  order  with  first  order  dependence  on 
Q:  and  Q,.  respectively. 

In  order  to  determine  the  order  of  reaction  and  thereby  the  rate  controlling 
step,  the  absorption  data  collected  at  2813  K under  various  partial  pressures  of  CO 
are  plotted  in  Fig.  5.14  for  carbon  and  Fig.  5.15  for  oxygen  using  the  integrated 
form  of  the  first  order  rate  equations  (5.50b)  and  (5.51b).  It  can  be  seen  that  in 
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Fig.  5.14  Carbon  absorption  kin 
equation 
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Fig.  5.15  Oxygen  absorption  ki 
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each  case,  ihc  data  are  split  into  four  distinct  lines,  each  representing  a different 
partial  pressure  of  CO.  Based  upon  the  same  reasoning  as  in  nitrogen,  this  behavior 
suggests  that  the  rate  controlling  step  in  C and  O absorption  in  liquid  Nb  is  not 
first-order,  implying  that  the  elementary  step  111  does  not  control  the  overall 

For  the  possibility  of  a second-order  reaction,  the  equation  (5.59)  and  (5.62) 
should  be  integrated  to  plug  the  data  into  it  directly.  Talcing  into  consideration  that 
the  ratio  of  O to  C is  constant  during  the  absorption  period,  both  equations  can  be 
transformed  into  the  equations  (5.63a)  and  (5.64a)  which  are  just  common  form  of 
second  order  rate  equations. 


dCc/dt  = k,'(A/V)(Cc,  Co,  - Cc  Co) 
= k,  (A/V)(  Cc’  - Cc*  ) 


(5.63a) 


or  in  integrated  form, 


lln  ( (£*Cc)  + ln  ( "ebre^1  = *■  4 1 (5  63b) 


dQ/dt  = k.'(A/V)  ( Co  Co.  - Cc  Co  ) 
= k.  (A/V)  ( C„,!  - Co1 ) 


(5.64a) 


integrated  form, 


(5.64b) 


where  k,  = 0.76  k,'  and  k,  = k»70.76. 

Thus,  the  integrated  forms  of  second  order 
(5.64b)  for  oxygen,  were  used  in  dealing 


equations,  (5.63b)  for  carbon  and 
i the  absorption  data.  Fig.  5.16  and 
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Fig.  5.16  Carbon  absorpiion  kinetics  correlated  wilh  second  order  rate 
equation 
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Fig.  5.17  show  the  second  order  plot  for  C and  O daia  ai  2*13  K in  which  il  can 
be  seen  dial  the  data  representing  the  different  partial  pressure  of  CO  fall  on  a 
single  straight  line  within  the  limits  of  experimental  error.  Thus,  il  can  be  said  that 
elementary  step  II  is  the  rate  controlling  step  for  C and  O dissolution  into  liquid  Nb 
from  CO  gas. 


controlled  by  dissociati 


second  order,  temperature  dependence  of  the  reaction  rate  can  be  derived  by 
examining  the  temperature  dependence  of  the  second  order  rate  constant,  k.  in 
equation(5.63b).  Fig.  5.18  shows  the  second  order  plot  for  the  data  obtained  at 
different  temperatures.  It  is  seen  that  the  data  collected  at  a given  temperature  falls 
on  a single  line  without  significant  deviation.  Taking  the  second  order  rate  constant 
at  each  temperaturc(from  the  slopes  in  Fig.  5.18),  an  Arrhenius  type  plot  of  those 
values  were  made  in  Fig  5.19.  Based  on  the  values  of  the  slope  and  the 
intercept(TabIe  10),  the  rate  equation  for  the  absorption  in  liquid  Nb  is  obtained  as 
(dCc/dtU  = 12.1(A/V)  exp(-26700/D(  Cc  C„  - Q C„  ) (5.65) 

where  Cc  and  Co  are  in  at  %,  A is  in  Cm',  V is  in  Cm’,  and  T is  in  K. 

It  is  to  be  noted  that  the  k.  value  obtained  in  this  plot  is  0.76  k.'  because  the  value 
of  O concentration  is  replaced  with  the  value  of  C concentration  in  calculation  for 
the  second  order  plot.  Since  O/C  ratio  remains  constant  at  0.76  during  the  entire 
dissolution  process,  the  oxygen  absorption  rate  can  be  represented  as 


(dCo/dtU  = 9,20(A/V)  exp(-26700/T)(  Cc  Q*  - Q C, ) ■ 


(5.66) 
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Fig.  5.18  Second  order  plot  of  carbon 


absorption  data 


Fig.  5.19  Arrhenius  type  plot  of  second  order  rate 
absorption  in  CO  atmosphere 
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Tabic  10.  Apparent  activation  energies  on  CO  -liquid  niobium  reaction 


slope  (Q/R)  intercept(ln  kj 


Absorption  26700  2.23 

Desorption  47500  10.55 


Fig.  5.19  and  Fig.  5.23. 
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The  Order,  of  ReaaiQn.agd  .Raie_Equation  for  CQ  Pcsorpiion 

In  the  analysis  of  desorption  data,  the  strategy  is  basically  analogous  to 
absorption  data  analysis  although  oxygen  desorption  docs  not  occur  through  oxide 
evaporation.  It  is  expected  that  desorption  rates  of  both  elements  are  the  same. 
However,  it  is  more  convenient  to  take  the  oxygen  data  instead  of  carbon  because 
desorption  process  through  CO  evolution  lead  to  carbon  residue  when  carbon  exists 
in  excess  at  the  beginning  of  desorption.  The  oxygen  desorption  data  in  Fig.  4.14  is 
plotted  in  Fig.  5.20  and  Fig.  5.21,  respectively,  in  accordance  with  the  integrated 
form  of  equations  (5.51b)  and  (5.64b).  The  first-order  plot  shows  a splitting  similar 
to  that  in  absorption,  while  the  data  are  correlated  well  with  the  second-order  plot. 
This  behavior  implies  that  the  desorption  reaction  is  again  controlled  by  elementary 

To  investigate  how  desorption  rate  is  varied  with  temperature,  the  second 
order  plot  is  drawn  in  Fig.  5.22  for  the  desorption  data  of  oxygen  obtained  at 
various  temperatures.  In  this  plot,  the  modified  integrated  form  of  second  order  rate 
equation,  equation  (5.67),  is  utilized  assuming  that  partial  pressure  of  oxygen  in 

1/Q,  - 1/Co,  = k(A/V)t  (5.67) 
Talcing  the  rate  constants,  k,  at  each  temperature  from  Fig.  5.22  (slopes),  the  value 
of  In  k is  plotted  against  IfT  in  Fig.  5.23.  Based  upon  the  values  of  slope  and 
intercept(Table  10),  the  rate  equation  for  the  desorption  process  can  be  expressed  as 
follows: 

(dCo/dt)*,  = (dQ/dt)_  « 3.8(W0,(A/V)  exp(-47500/T)(  Q,  Co.  - Q CD  ) (5.68) 


Cc  and  Co 
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Fig.  5.20  First  onler  plot  of  oxygen  desorption  kinedcs  in  atmospheres 
with  various  partial  pressures  of  CO 
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Second  order  ploi  of  oxygen  desorption  kinetics  in  atmospheres 
with  various  partial  pressures  of  CO 


Fig.  5.21 
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Fig.  5.22  Second  order  plot  of  oxygen  desorption  data  in  pure  Ar  atmosphere 
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Fig.  5.23.  Plot  of  Ln  k,  second  order  rate  constant  of  oxygen  desorption  in 
pure  Ar,  versus  1/T 


CHAPTER  6 
CONCLUSION 

Using  an  electromagnetic  levitation  technique,  absorption  and  desorption 
behaviors  of  nitrogen  and  carbon  monoxide  in  liquid  niobium  were  investigated 
under  various  experimental  conditions  of  temperature,  partial  pressure,  gas  (low  rate, 
and  the  processing  time.  Some  important  results  are  as  follows: 
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respectively,  implying  that  ihe  overall  reaction  rale  is  controlled  by  ihc 
subslep  of  dissociative  adsorption  (or  recombination  of  C and  O in  the  case 
of  desorption). 

The  rate  equations  for  C/O  absorption  and  desorption  are  found  to  be 
(dCVdt).  = (l/0.76)(dCVdtU  = lll(A/V)exp<-267Wr)(Cc.  Co.  - Q Co) 


(dCVdt)»  = (dCo/dt)u  = 3.8<M0*(A/V)  cxp(-47500/T)(  Q.  Cc  - Q Co  ) 
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